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PREFACE 

This  review  Is  an  attempt  to  summarize  the  present  sUte  of  knowIed|;e  of  the 
process  by  which  a  seismic  wave  is  ^nerated  by  an  underground  explusioii.  An  at¬ 
tempt  has  been  made  to  start  with  the  highly  compressed  volume  of  gas  formed  at  the 
completion  of  detonation  and  proceed  step  by  step  to  the  seismic  signal.  Research 
investigators  in  this  field  will  recognize  immediately  that  there  are  gaps  in  our  knowl¬ 
edge  about  critical  parts  of  this  process  that  m»st  be  bridged  by  reasonable  assump  ¬ 
tions  about  material  behavior. 

The  author  is  a  selsmul'vglst,  and  this  report  is  directed  primarily  to  other  seis¬ 
mologists.  For  this  reason,  experts  in  shock-wave  physics  may  find  that  the  sections 
on  the  highly  stressed  region  around  the  explosion  contain  too  much  detailed  discussion 
of  fundamental  principles,  without  a  corresponding  coverage  in  the  sections  on  elastic 
waves. 

The  author  wishes  to  acknowledge  with  gratitude  the  invaluable  assistance  given 
him  in  the  preparation  of  this  report  by  the  staff  of  VESIAC,  Institute  of  Science  and 
Technology,  The  University  of  Michigan.  In  particular,  Mr.  Thomas  Caless,  Mr. 
Robert  Haven,  and  Mrs.  Elaine  Medor  have  provided  services  without  which  the  com¬ 
pletion  of  this  work  would  have  been  ver.  difficult. 

Finally,  the  author  wishes  to  express  his  gratitude  to  his  colleague.  Prof.  Otto 
W.  Nuttli,  Department  of  Geophysics  and  Geophysical  Engineering,  St.  Louis  Uni¬ 
versity,  who  read  the  manuscript  during  its  preparation  and  offered  many  helpful 
suggestions. 
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THE  GENERATION  OF  THE  PRIMARY  SEISMIC  SIGNAL 
BY  A  CONTAINED  EXPLOSION 

ABSTRACT 

In  order  to  analyze  the  mechanism  by  which  buried  explosions  generate  >eismic 
waves,  the  region  around  the  detonation  is  divided  into  three  zones:  (1)  a  strong- 
shock  (hydrodynamic)  zone,  (2)  a  transitional,  nonlinear  zone,  and  (3)  the  elastic 
region.  Experimental  equatlon-of-state  data  are  used  in  calculating  the  history  of 
the  propagating  stress  wave,  because  there  is  as  yet  no  comoiete  theory  of  the  re¬ 
sponse  of  solids  to  rapidly  applied  stresses  which  exceed  their  yield  and  crushing 
strengths.  Results  indicate  that  the  proper  equations  of  motion  are  known,  as  rre 
analytic  and  numerical  methods  for  solving  them.  The  theory  of  shock  waves  in  a 
fluid  is  applicable  to  the  close-in  region,  and  srvetal  failure  mec.hanlsms  have  been 
postulated  for  the  transition  region. 

The  objective  of  the  analysis  is  the  determination  of  the  stress  waveform  at  the 
inner  boundary  of  elastic  behavior.  The  peak  amplitude  and  spectral  content  depend 
on  the  yield,  the  type  of  medium,  and  the  ambient  stress  (depth  of  burial).  Tiie  ex¬ 
perimental  determination  of  the  shortest  range  at  which  elastic  behavior  begins  is 
difficult  in  principle,  because  solutions  valid  at  great  distances  are  not  applicable. 

Experimental  determinations  of  the  amplitude-yield  relationship  must  take  into 
account  (1)  the  combined  effect  of  the  shift  of  the  spectral  peak  to  lower  frequencies 
as  the  yield  increases  and  (2)  the  low-pass  filtering  properties  of  the  earth.  Em¬ 
pirical  studies  of  the  effect  of  the  shot- point  medium  are  influenced  by  the  properties 
of  the  explosive  material.  The  effect  of  source  depth  is  difficult  to  Isolate  because 
in  most  test  sites  the  medium  changes  with  depth. 


I 

INTRODUCTION 

Of  the  three  physical  processes  involved  In  seismic  ex¬ 
ploration,  namely  the  initiation  of  the  seismic  waves,  their 
propagation,  reflection,  refraction,  and  dispersion,  and  the 
recording  of  some  function  of  the  motion  of  the  surface,  we 
possess  the  least  satisfactory  understanding  of  the  initiation 
process. 

J.  A.  Sharpe,  1942 

An  assumption  underlies  the  entire  research  program  on  the  Identification  of  underground 
nuclear  explosions  by  the  seismic  signal  recorded  at  suitably  placed  stations;  that  Is,  such  ex¬ 
plosions  represent  sources  having  unique  properties  which  arc  adequate  to  dlstingiilsh  them  from 
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all  other  seismic  sources.  The  discovery  of  these  properties  will  permi*  the  establishment  of 
criteria  for  deciding  that  a  particular  seismogram  oi  set  of  EclsmoRrams  could  have  been  pro¬ 
duced  only  by  an  explosion. 

The  theory  of  wave  propagalloi'.  In  the  earth  has  progressed  to  the  point  where  It  would  be 
possible  to  compute  the  seismogram  resulting  from  an  Input  of  kt'.jwn  form  at  the  source,  given 
the  detailed  properties  of  the  geologic  structure  between  a  seismic  source  and  observing  point, 
and  the  response  characteristics  of  the  seismograph.  Theoretical  seismograms  for  simplffind 
structures  are  being  computed  for  a  number  of  applications.  Of  course,  not  enough  is  known  of 
the  constitution  of  the  cn^st  cr  mantle  to  make  this  ;ossible  at  present,  but  ambitious  research 
efforts  now  in  progress  promise  to  yield  the  required  data  eventually. 

An  excellent  example  of  current  ability  to  derive  a  particular  wave  on  the  seismogiam  from 
a  detonation  of  known  yield  is  offered  by  the  rK>mputation  of  arrivals  for  tnree  nuclear  explo¬ 
sions  by  Werth  et  al.  [1].  Their  methods  and  results  indicate  both  the  potential  of  the  analytic 
tools  available  and  the  severe  limitations  of  the  data  available  to  perform  the  computations. 

The  progress  that  has  been  made  in  the  computation  of  theoretical  seismograms  makes  the 
question  of  the  properties  of  the  source  a  pressing  one,  for  we  are  asking,  "What  should  be  the 
input  to  our  computation?" 

All  efforts  to  answer  this  question  have  begun  by  surrounding  the  source  with  a  sphere  with 
a  radius  targe  enough  that  all  processes  outside  of  this  sphere  can  be  described  by  infinitesimal 
strain  theory,  either  elastic  or  allowing  for  frictional  losses.  Then  the  question  becomes  one  of 
the  proper  loading  (or  set  of  stresses)  to  distribute  over  this  surface  to  represent  the  output  of 
the  explosion  in  the  form  ot  seismic  energy. 

This  report  will  summarize  present  knowledge  of  the  seismic  loading  from  a  contained  ex¬ 
plosion.  Ofparticular  interest  are  the  effects  on  the  resulting  signal  of  the  yield  of  the  explosive  de¬ 
vice,  the  medium  in  which  it  is  buried,  and  the  depth  below  the  earth's  surface.  Because  small- 
and  large-scale  geologic  factors,  which  are  not  an  intrinsic  part  of  the  wave-generating  process,  in¬ 
fluence  the  character  of  the  signal  profoundly  even  at  moderate  distances,  only  the  close-in  effects 
will  be  considered — that  is,  effects  at  the  shortest  ranges  at  which  the  signal  is  carried  by  elas¬ 
tic  waves.  For  example,  Adams  et  al.  [2]  have  found  (at  the  Nevada  Test  Site)  that  inhomogeneitles 
in  the  medium  seriously  disturb  measurements  of  ground  motion  at  scaled  ranges  beyond  800  to 
1000  feet.' 

'Scaled  range  »  range/fW(kllotona)]  ;  see  the  appendix  for  dlscusslnn  of  scaling. 
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A  great  deal  of  Important  work  on  close-in  effects  of  explosive  sources  Kis  been  done  by 
researchers  who  are  Interested  primarily  in  weapons  eitects  and  protectb'.  construction.  For 
obvious  reasons,  much  of  this  work  has  dealt  with  surface  or  air  bursts,  because  (^e  loading 
due  to  the  Intense  air  blast  constitutes  an  Important  part  of  the  source  mechanism  for  these  shots, 
the  results  of  this  work  will  not  be  considered  here.  Only  contilneo  explosions  are  of  Interest 
here,  and  throughout  this  paper  the  term  "explosion"  will  be  synonymous  with  "contained  e/:* 
plosion"  unless  otherwise  Indicated.  Minor  venting  often  occurs  when  burlea  exploutons  are  de¬ 
tonated;  an  explosion  will  be  considered  contained  as  long  as  the  associated  air  blast  has  a  negli¬ 
gible  effect  on  ground  motion.  Any  treatment  of  viscous  or  other  losses  in  the  seismic  waves 
during  propagation  is  also  outside  the  scope  of  this  report.  There  are  imporUnt,  as  representing 
a  significant  effect  of  the  medium  through  which  the  signal  travels.  This  discussion,  howc.  ur, 
will  concern  only  those  loss  mechani»m.s  which  represent  an  effect  of  the  medium  in  which  the 
explosion  occurs. 

Until  very  recently,  most  of  the  work  on  the  factors  affecting  the  seismic  output  of  an  ex¬ 
plosion  was  done  by  persons  who  were  interested  in  the  use  of  explosives  as  controlled  energy 
sources;  exploration  seismologists,  seismologists  working  with  quarry-  and  mine -blast  oper¬ 
ations,  and,  to  a  lesser  extent,  demolition  experts. 

Seismic  exploration  detonation  conditions  pertain  more  closely  to  the  underground  nuclear 
test  Identification  and  detection  problem  than  do  quarry  blasts.  Even  though  the  size  of  a  large 
quarry  or  construction  bias*  is  somewhat  closer  to  the  scale  of  a  nuclear  explosion,  the  fact  that 
the  shot  is  always  near  a  free  face  and  is  desijy^d  to  move  rock  means  that  a  quarry  shot  is 
never  contained.  This  fact  significantly  alters  t. »  seismic  output.  Furthermore,  large  commer¬ 
cial  blasts  are  always  distributed  in  space  and  time  through  the  delayed  firing  of  a  number  of 
holes.  For  these  reasons,  the  results  of  theoretical  and  experimental  research  carried  out  in 
connection  with  the  use  of  explosives  in  exploration  are  of  more  value  for  this  study  than  those 
based  on  observations  of  quarry  blasts.  However,  some  excellent  basic  work  on  the  mechanism 
of  rock  breakage,  especially  that  carried  out  by  the  U.  S.  Bureau  of  Mines,  bears  directly  on  the 
subject  of  this  report. 

Because  of  the  differences  between  chemical  and  nuclear  detonations,  research  based  on  the 
study  of  high  explosives  alone  cannot  possibly  solve  this  problem.  Consequently,  only  very  recent 
or  In-progress  research  (theoretical  and  experimental)  on  nuclear  sources  will  aid  in  elucidating 
the  wave-generating  processes.  Most  of  the  Important  information  on  the  behavior  of  earth  ma¬ 
terials  in  the  highly  stressed  region  around  the  explosion  Is  still  to  be  gathered.  As  this  report 
was  being  written,  an  extensive  series  of  underground  explosions  Is  In  progress  under  the  spon¬ 
sorship  of  the  U.  S,  Atomic  Energy  Commission;  It  is  very  likely  that  analysis  of  data  from  these 
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experiments  will  provide  numerical  values  to  be  used  In  the  theoretical  results*  outlined  here,  at 
least  for  the  selected  materials  at  the  test  sites. 

2 

THE  MODEL  OF  AH  EXPLOSION  AS  A  SEISMIC  SOURCE 

If  it  were  possible  to  surround  an  explosion  with  a  surtkce  of  known  shape  on  which  the  com¬ 
plete  stress  history  at  each  point  was  known  as  a  function  of  the  parameters  describing  the  ex¬ 
plosion  (yield,  type  of  explosive,  medium  around  the  charge,  depth  of  burial),  and  if  this  surface 
were  so  chosen  that  only  elastic  processes  took  place  outside  of  it,  the  differential  equations  de¬ 
scribing  elastic  wave  propagation  could  be  integrated  to  give  the  ground  motion  at  any  point  in  the 
neighborhood  outside  this  source  region. 

The  theory  of  wave  generation  by  explosions  has  advanced  through  .  succession  of  efforts  to 
get  an  Increasingly  accurate  picture  of  the  loading  function.  The  earliest  efforts  [3,  4,  5]  simply 
ignored  the  complex  processes  Inside  the  hypothetical  surface,  and  proceeded  from  an  Intelligent 
conjecture  of  the  nature  of  the  loading.  This  approach  made  it  impossible  to  relate  actual  ampli¬ 
tudes  and  waveforms  of  ground  motion,  stress,  or  strain  to  the  yield,  but  made  possible  a  qual¬ 
itative  picture  of  the  manner  in  which  the  explosion  parameters  a<i'ect  the  results  to  be  obtained. 

The  next  stage  was  to  discuss  qualitatively  the  effects  of  the  nonlinear,  irreversible  region 
[6,  7,  8,  9] .  Some  early  studies  have  been  followed  up  by  the  present  concerted  efforts  to  treat 
in  detail  the  processes  occurring  in  the  highly  stressed  region.  Authorities  generally  agree  on 
how  to  handle  the  "strong"  shock  zone  immediatoiy  around  a  nuclear  detonation;  and,  as  mentioned 
above,  there  is  no  serious  problem  in  the  elastic  'cne.  It  is  true  that  loss  mechanisms  at  in¬ 
finitesimal  strains  must  be  taken  into  account  If  the  seismic  waveform  is  to  be  accurately  pre¬ 
dicted  [10].  The  greatest  uncertainty  lies  in  the  transition  zone,  where  stresses  are  falling  to 
the  neighborhood  of  the  crushing  strength  and  then  the  yield  strength  of  the  medium.  However,  as 
pointed  out  by  Fugelso  [11,  p.  1],  the  unknown  factors  in  the  identification  of  a  seismic  event  as 
a  blast  lie  in  this  region.  The  various  approaches  outlined  in  this  report  can  almost  be  categorized 
by  the  number  of  subregions  into  which  this  zone  of  transition  is  divided. 

Cushing  and  Relly  [12]  suggest  a  useful  coarse  division  of  the  zones  around  the  shot.  They 
refer  to  the  close-in  nonlinear  region,  in  which  shock  pressure  far  exceeds  the  crushing  strength 
and  shear  stresses  are  negligible  compared  to  compressive  stresses,  and  to  the  terminal  non¬ 
linear  region,  in  which  the  shock  pressures  are  at  or  above  the  elastic  limit  of  the  medium,  and 
are  generally  comparable  to  the  medium's  crushing  stress,  and  in  which  the  shear  strength  must 
be  considered.  It  is  possible  and  desirable  to  define  the  limits  of  these  regions  quantitatively,  and 
to  further  subdivide  the  terminal  nonlinear  region  [13], 
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This  report  will  present  the  theory  of  wave  generation,  starting  from  the  Immediate  neighbor¬ 
hood  of  the  expluslon  and  proceeding  outward,  Th?  main  features  of  the  ci.^c-n-in  nonlinear  region 
will  be  discussed  first.  Then  the  various  vlewi>ofnt8  on  the  transition  region  will  !,•>  reviowe'l, 
and  finally  the  character  of  the  seismic  signal  at  the  beginning  of  elastic-wave  propagation  will  be 
described.  Present  knowledge  of  the  waveform  (time  history)  or  spectrum  (frequency  content)  of 
the  input  seismic  signal  will  be  summarized.  The  report  tri'l  also  review  field  data  from  chem¬ 
ical  and  nuclear  explosions,  and  will  present  the  various  empirical  scaling  rules  t;at  ha^e  been 
developed,  compared  with  the  theoretical  results. 

The  physics  of  the  explosion  itself  will  not  be  discussed.  Detailed  discussion  of  the  detona¬ 
tion  process  for  chemical  explosions  may  be  found  in  Cole  [14],  Cook  [15],  and  Jones  and  Miller 
[16] .  It  will  be  assumed  that  the  detonation  results  in  an  initial  spherically  symmrlrlc  mtense 
shock  front.  Evidence  of  departures  from  symmetry  of  the  seismic  .  ignal  from  an  explosion  will 
be  discussed  in  Section  6.4. 


3 

SHOCK  WAVS  PROPAGATION  AND  THf  HUGONIOT  EQUATION  OF  STATE 

3.1.  THE  RANKINE-HUGONIOT  EQUATIONS 

Stresses  in  the  immediate  vicinity  of  the  explosion  reach  levels  of  hundreds  of  kilobars  or 
megabars,  and  exceed  the  strength  of  any  earth  material  enough  that  the  rigidity  is  completely 
negligible  and  the  phenomena  may  be  described  by  treating  the  medium  as  a  perfect  fluid.  Studies 
that  attempt  to  work  in  detail  from  the  explosion  outward  treat  the  shock  wave  as  propagating  in 
a  fluid.  The  phenomena  in  this  case  are  governed  by  the  Rariklne-Hugcniot  equations  for  a  shock 
wave  in  a  fluid.  These  equations  will  be  reviewt\l  here  because  they  are  basic  to  much  that  fol¬ 
lows.  More  detailed  treatments  of  these  relations  may  be  found  in  worksdevoted  to  the  present 
problem  [13,  18,  20]. 

The  Rankine-Hugoniot  relationships  are  derived  from  the  principles  of  conservation  of  mass, 
momentum,  and  energy  across  a  steep  front  at  which  pressure,  density,  and  particle  velocity  are 
discontinuous.  If  the  shock  front  is  advancing  with  a  speed  U  into  an  undisturbed  region  of  fluid 
in  which  the  density  is  and  the  ambient  pressure  p^,  conservation  of  mass  across  the  front 
requires  that  the  particle  velocity  u  and  the  density  p^  Immediately  behind  the  front  satisfy  the 
equation  p^(U  -  u )  ■  p^U,  which  may  be  written  as 


where  v  •  1/p  is  the  specific  volume. 
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Conservation  of  momentum  is  used  to  relate  the  pressure  Jump  across  the  front  to  the  other 
2  ? 

parameters:  p  -  p  «  p  U  -  p  (U  -  u)  .  If  p  is  eliminated  by  using  Equation  1  and  calling  the 

B  iJ  O  8  B 

overpressure  p^  -  p^  ■  P(  then 

p=p^Uu  (2) 

Conservation  of  energy  requires  that  the  work  done  by  t:;e  pressures  at  the  front  equals  the 
change  In  the  total  energy  of  the  fluid.  Referred  to  a  unit  mass  where  represents  the  Initial 
Internal  energy  and  the  change  at  the  front,  the  total  change  In  energy  per  unit  mass  U>i  • 

ternal  and  kinetic),  assuming  an  adiabatic  process,  is  (1^  ”  5  Equating  this 


to  work  per  unit  mass. 


p^U.p^(U.u) 


,  yields 


Po  o  5 


77  +  1  ^  5  (U  -  u)^.  This  can  be 
p^  s  2 


rewritten,  using  Equation  2: 


-P8(^-^)"Ps^'^0-  V 

These  three  basic  relations  can  be  combined  in  various  ways  to  obtain  the  desired  shock 
parameters  In  terms  of  those  that  are  more  easily  measured.  The  velocity  of  shock  propagation 
and  the  particle  velocity  In  terms  of  pressure  and  specific  volume  are  found  from  Equations  1 
and  2: 


2  Pb  *  P.' 

V..  -  V.  o 

O  8 


]  “^“^Ps'Po'^^,- V 

From  Equations  5  and  3a,  the  change  in  internal  energy  Is 

<*s-V*I<Ps*Po^^'^o- V 

Equation  6  is  called  the  Hugonlot  equation  of  the  fluid  [IT].  It  is  the  locus  of  all  final  states  p^, 
v^  that  can  be  attained  from  an  Initial  state  p^,  v^  by  passing  a  shock  through  the  fluid. 

The  significance  of  the  various  expressions  derived  above  can  be  visualized  through  the 
pressure-volume  diagram  In  Figure  1  (see  Reference  13,  Figure  4.1,  Reference  18,  Figure  1, 
and  Reference  19,  Figure  3.3).  From  Equation  3a,  the  total  energy  per  unit  mass  behind  the 
shock  In  p  (v  -  v  ),  the  area  of  the  rectangle  bounded  by  p  ■  p  ,  v  ■  v  ,  and  v  ■  v  .  From 
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FIGURE  1.  ENERGY  RELATIONS*  ON  THE  GRAPH  OF  THE  EQUATION 
OF  STATE  (After  BIsho;-,  Cf.aBzeyka,  and  PorzeT) 

Equation  5,  the  kinetic  energy  per  unit  mass  is  5  (p  -  p  )  (v  -  v  ),  the  area  of  the  triangle  above 

B  O  O  B  j 

the  line  of  shock  compression.  From  Equation  6,  the  area  of  the  trapezoid  ^  (p^  +  p^)  (v^  ~  v^) 
represents  the  Internal  energy.  The  area  under  the  expansion  curve  back  to  the  initial  state 
represents  fpdv  work,  and  is  available  as  potential  or  hydrodynamic  energy  which  is  given  up 
during  expansion.  The  area  between  the  compression  line  and  the  expansion  curve  represents 
energy  loot  from  the  shock  by  the  increase  in  entropy  across  the  steep  shock  front.  Because  it 
is  no  longer  available  to  support  the  propagation  of  the  shock,  this  energy  has  been  called  waste 
heat.  It  appears  as  a  final  rise  in  the  temperature  of  the  medium  [19].  Heat  conduction  occurs 
so  slowly  that  the  deposition  of  waste  heat  is  permanent  as  far  as  the  shock  wave  is  concerned. 
For  solids,  and  less  exactly  for  fluids,  the  coefficient  of  thermal  expansion  is  small  enough  ttuit 
the  final  volume  in  close  to  v^  [13].  The  equation  of  state  or  Hugonlot  curve  Is,  of  course,  a 
property  of  the  material  and  at  present  is  Iwisod  on  empirical  values  of  p  vs  v. 
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3.2.  THE  VARIATION  OF  U  WITH  p 

A  useful  property  of  shock  wave  propagation  can  be  derived  from  Equation  4  by  differentiating 
with  respect  to  p^: 


1 


dU 

w/s: 


\v  -  V  ; dv 
_ \  o  8/  s 

(v  -  V  (p  -  p  ) 

O  S  ^8  O 


‘*‘*8 

The  quantity  i8  obtained  from  the  Hugoniot  eq^'*tion  of  8tate  for  the  medium  and  is  a  function 
of  V  .  The  shock  velocity  U  will  Increase  with  increasing  pressure  as  long  as 


(7) 


For  all  materials,  both  sides  of  the  Inequality  as  written  are  positive  (pressure  decreases  with 
Increasing  specific  volume),  and  the  Inequality  Is  satisfied  If  the  line  Joining  the  shocked  state  to 
the  Initial  state  makes  a  smaller  angle  with  the  negative  v-axis  than  does  the  tangent  to  the  equa* 
tlon  of  state  at  the  point  (p  ,  v  ).  As  shown  by  Rice  et  al.  fl7]  and  discussed  by  Cushing  and  Reily 
[12],  if  the  equation  of  state  has  a  cusp  (as  It  will  in  any  solid  with  a  non-zero  yield  stress),  and 
a  somewhat  greater  crushing  stress,  there  will  be  a  region  In  which  shock  velocity  increases 
with  decreasing  pressure,  and  in  place  of  a  stable  single  shock  front,  a  two-wave  structure  forms. 
This  will  occur  when  the  shock-front  pressure  haa  decayed  to  the  neighborhood  of  the  crushing 
stress  of  the  medium  In  the  terminal  nonlinear  retb>n.  The  shock  front  velocity  can  become  less 
than  the  velocity  of  the  elastlu  compressional  wave,  so  that  an  elastic  precursor  runs  ahead  of 
the  shock  front. 


8.8.  INTRODUCTION  OF  ADIABATIC  LAW  [13,  20] 

The  internal  energy  can  be  eliminated  from  the  equation  of  state  (Equation  6)  by  Introducing 
the  law  for  adiabatic  expansion  of  an  ideal  gas.  The  change  in  Internal  energy  of  the  system  in 

f''f 

going  from  an  Initial  to  a  final  state  Is  I^  -  I^  >  I  pdv,  where  the  integration  is  taken  along  an 

Jv, 

adlabat,  pv^  ••  const.;  >  *  ratio  of  specific  heats.  Integrating, 


»f*^ 


Vf- Vi 

y  -  1 


(7a) 
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Bishop  takes  the  value  of  for  the  adiabat  connectinfr  the  ambient  and  shocked  states 

and  equates  it,  as  an  approximation,  to  the  right  hand  side  of  Equation  6,  obtaining 


(y  •  t)"^  (p  V  -  p  V  )  -  i  (p  +  p  )(v  -  V  ) 
a  *^8  8  *^o  o  2  *^8  o  o  8 


(8) 


He  shows  that  this  approximation  holds  for  strong  shock.'*  In  which  p  »  p  by  obtaining  for  this 

s  o 

p_  y-  + 1 

case —  - - j,  a  result  given  by  Taylor  [21],  He  also  shows  that  Equation  8  is  vail-  for  very 

^o  ’'s  '  ^ 


low  pressures  (acoustic  case),  by  demonstrating  that  y_  which  is  the  differential  form 

of  the  equation  of  the  adiabat. 


s 


Therefore,  there  is  justificatlmi  for  proceeding  to  use  an  ’’effectlve-a/erage"  adiabatic  ex¬ 
ponent,  y  ,  in  any  "hydrodynamic”  medium.  Other  authois  have  proccoJed  to  use  this  y  without 
explicitly  examining  the  approximation  involved  in  Equation  8.  By  replacing  y^,  the  exponent  at 
p  ,  with  y,  a  value  for  states  in  the  neighborhood  of  p  ,  and  using  the  Ranklne-Hugontot  relations, 
several  results  given  by  Taylor  are  readily  obtained. 

First,  from  Equation  8, 


„  (y-t)r(yrl)-2- 

fs___ _ ^ 

Po  Ps 

°  (yr  1)  r  (y.  1).^ 

Po 


(9) 


Introducing  c 


2 


(the  acoustic  velocity  in  the  undisturbed  medium)  and  Equation  9  into 


Equations  4  and  5  gives 


(10) 


U 

8 

IT 


y.  1  +(y+  !)■ 


(11) 


Bishop  considers  a  medium  to  l>e  in  a  "hydrodynamic"  state  if  y  is  a  slowly  varying  function 
of  pressure. 
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f 

I 

t 

'I 

i 


! 


A  shock  Is  considered  "strong"  If  Is  sufficiently  greater  than  so  that  Equations  9  through 
11  take  the  form  \2l] : 


P3  y  1 

(9a) 

y  ♦  1  Ps 

~f“~2rr 

c  *^o 

(10a) 

0 

“s  2 

(11a) 

A  very  useful  form  of  the  equation  of  state  Is  obtained  by  Introducing  ®  (f’le  clia*.ge 

In  specific  volume)  and  P  ■  Pg  •  Pq  overpressure)  Into  Equation  9. 


o  Ps  Ps  -  1 

°  (y  -  1)  V  (y  ♦  1)  ^  2y  +  (y  +  1) 


2y  +  (tv* 

*^o 


or 


(12) 


Bishop  interprets  yp^  as  the  adiabatic  elastic  modulus  for  a  condensed  medium,  where  y  is 

yPo  2 

at  the  ambient  state.  If  o  Is  the  ambient  compresslonal  wave  velocity,  then - ■  c  ,  so  that 

°  Po  ° 

2 

yp^  ■  p^c^  ■  K^,  an  elastic  modulus.  Therefore,  Equation  12  becomes 


>2  (y  +  1)  +-;r 


P 


(13) 


-  • 
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K  p/K 

Dividing  through  by  ,  --y~  =  1  +  1/i  (y  +  1) 


(14) 


If  the  overpressure  p  Is  small  compared  to  the  elastic  constant  K^,  Vj^  is  krU  Equation 

13  becomes 


»  1,  or  K 


P  ^  P 

b  V. /v  -Av  /v 

n  o 


This  merely  confirms  the  Interpretation  of  as  an  elastic  modulus,  i.e.,  the  ratio  of  small 
stress  to  small  strain.  Then,  for  large  stress.uS, 


--Pr-  =  K  (1  *  X)  « 

•'h'  •'o  ®  ® 


where 

X  -  1/2  (y  ♦  1)  jH.  -  (Kg/K^)*^ 
o 


(15) 


(16) 


and  K  is  the  elastic  modulus  for  finite  strains.  The  quantity  x  is  a  factor  to  correct  the  Infin- 
s 

Itesimal  strain  modulus  to  a  finite  strain  modulus,  and  is  a  function  of  the  stress. 


4 

In  notation  more  familiar  to  seismologists,  c^  is  the  P-wave  velocity,  K=k+-5M»X+2/i, 

O  O  o 

in  which  k  is  the  bulk  modulus,  p  the  rigidity  modulus,  and  A  is  Lame^'s  constant.  If  the  velocity 
of  propagation  of  elastic  shear  waves  is  designated  by  c^,  then 


1  - 


liL  =  i  .1 
3  K  3 

o 


From  Equations  1  and  2,  Equation  13  becomes 


1/2  (y  +  1) Cjj*  (uU)*^  (17) 

y  can  then  be  determined  from  the  experimental  propagation  data,  U,  u,  and  c^.  From  Equations 

U2 

16  and  17,  another  form  for  x  Is  x  »  — j  -  1. 

*^o 

Following  Bishop,  we  shall  use  the  value  of  x  to  describe  the  regime  in  which  the  processes 
are  taking  place;  I.e.,  when  x  Is  very  large,  the  shock  is  strong,  and  when  x  Is  very  small,  the 
behavior  Is  elastic.  Thus  transitional  states  can  bo  specified  for  a  given  medium  by  values  of  x. 
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If  observational  data  are  available.  If  overpressure  Is  used  to  specify  tbo  regime,  Equation  12 
can  be  expressed  as 


y  +  1 


'I 

l| 


2yx 
Y  *  1 


(18) 


Bishop  has  assembled  values  for  y,  based  on  Equation  17,  ioi  v.^rlous  media  from  data  col 
lected  by  several  Investigators,  These  values  are  presented  in  Table  1. 


For  X  greater  than  1,  y  is  generally  found  to  be  a  slowly  var]ding  function  of  pand,  from  Equa¬ 
tion  18,  overpressure  varies  linearly  with  x. 


The  equations  of  state  for  tuff  and  halite,  taken  from  Bishop,  are  shown  In  Figures  2  and  3. 
These  summarize  his  analysis,  and  will  discussed  in  detail  below.  On  these  figures  p^  ■  p^, 
the  ambient  pressure. 


TABLE  I.  ADIABATIC  EXPONENTS  FOR  SOME  GEOLOGIC 
MATERIALS  AND  FOR  WATER  [13] 


Material 

Density 

Adiabatic 
Exponent  Y 

Acoustic 
Velocity  c^ 

Granite  and  Basalt 

2.67  (gm/cm^) 

1.7 

4.9  (km/sec) 

Halite 

2.15 

2.1 

4.40 

Wet  tuff 

1.85 

3.2 

2.44 

Dry  tuff 

1.7 

2.1 

2.44 

Alluvium 

1.6 

2.7 

1.22 

Water 

p  =  30  kb 

5 

1.5 

Water 

p  -  100  kb 

3.8 

1.5 

Water 

p  ■  400  kb 

2.8 

1.5 
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WAVE  FRONT  OVERPRESSURE  (bars) 

FIGURE  2.  EQUATION  OF  STATE  FOR  HALITE  (p^s  -  60  bars)  (After  Bishop) 
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4 

THC  iTROMG  SHOCK  REGION 

According  to  the  qualitative  classification  of  Cushing  and  Relly  [12],  the  clos»-in  nonlinear 
region  Is  that  In  which  a  hydrodynamic  shock-wave  analysis  Is  valid.  A»  slated  above,  this  re¬ 
quires  that  shock  pressures  greatly  exceed  crushing  pressure  a.;i  that  shear  stresses  be  neg¬ 
ligible.  The  formidable  task  of  determining  the  thermou^namlc  properties  of  the  strongly 
shocked  material  must  then  be  accomplished.  The  objective,  from  the  viev’iJolni  of  this  report, 
is  to  determine  where  the  transition  to  the  terminal  nonlinear  region  occurs  and  how  much 
energy  remains  in  the  shock  wave  at  that  radius  f22,  p.  75] . 

The  Ranklne-Hugonlot  equation  of  state  is  required  for  the  medium;  at  present,  this  infor¬ 
mation  is  primarily  empirical.  The  equation-of- state  work  on  metals  gives  useful  tnslght  Into 
the  problem,  but  little  has  been  done-  on  porous  materials.  The  work  of  Chaszeyku  [22]  and 
Chaszeyka  and  Porzel  [18]  seems  to  be  the  best  available  on  materials  resembling  earth  ma¬ 
terials;  work  is  in  progress,  however,  by  Fugelso  [11]  and  Bishop  [13], 

Fugelso  points  out  that,  owing  to  the  characteristics  of  earth  materials  (l.e.,  inhomogeneity, 
complex  mineraloglcal  structure,  porosity,  and  small-scale  structural  defects),  there  is  no 
hope  of  finding  an  analytic  equation  of  state  to  fit  the  experimental  data  over  the  entire  range 
of  interest. 

Bishop  defines  the  hydrodynamic  region  as  that  In  which  the  effective  adiabatic  exponent  y 
varies  slowly  with  pressure.  The  lower  limit  of  the  hydrodynamic  range  is  a  critical  pressure 
at  which  crushing  ends  and  hydrodynamic  flow  begins.  Because  this  pressure  falls  in  a  region 
of  plastic  flow,  it  is  somewhat  indeterminate.  Bisi^op  approximates  this  critical  value  p^^,  the 
hydrodynamic  termination  pressure,  with  a  value  Pg  which  corresponds  to  a  value  of  x  ■  Xg  ■ 

0.5.  For  p  >  p  ,  the  x  -  p  curve  (see  Figures  2  and  3)  is  almost  a  straight  line  (y  -  constant), 
except  for  values  close  to  p^,  where  a  transition  zone  of  large  plastic  deformations  occurs. 

A  shock  Is  then  defined  as  being  "strong"  initially  If  x  >  100  for  any  radial  distance.  This 
occurs  for  nuclear  explosions,  but  never  for  chemical  detonations.  Figure  3  shows  that  x  ■  100 
for  tuff  at  about  1.5  megabars. 

Porzel  [19]  maintains  that  the  waste  heat  is  the  controlling  mechanism  for  a  well  contained 
spherically  symmetric  explosion,  such  as  the  RAINIER  event.  This  viewpoint  is  supported  by 
analysis  of  soll-Iike  mixtures  of  air  and  solid  particles  [18,  22],  Waste  heat  was  defined  above 
as  heat  which  remains  in  the  material  after  it  expands  adiabatlcally  back  to  the  ambient  pres- 
sure  following  the  passage  of  the  shock  wave,  Chaszeyka  has  shown  [22]  that  for  even  small 
amounts  of  air  Included  in  the  soir'or  rock-like  material,  virtually  all  of  the  energy  of  shock 
compression  goes  Into  waste  heat,'^ 
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For  all  materials,  there  is  some  pressure  value  high  enough  to  he?t  and  vaporize  the  ma¬ 
terial  irreversibly,  and  leave  It  In  a  gaseous  state  even  after  expansion  to  the  ambient  pressure. 
All  of  the  material  Inside  the  shock  radius  corresponoitig  to  this  "vaporization  pressure"  be¬ 
comes  part  of  the  cavity.  For  very  strong  shocks,  there  will  be  energy  In  the  v.’aste-hcat  part 
of  the  p-v  diagram  in  excess  of  that  required  to  vaporize  the  materlr.i-  this  will  appear  as  super¬ 
heat  In  the  gas.  Since  this  energy  Is  available  by  expansion  -1  the  gas,  it  Is  not  part  of  the 
waste  heat  defined  above.  Thus,  for  pressures  far  In  excess  of  the  vaporizallon  prcss'ire.  all 
materials  are  good  transmitters  of  shock  energy  fl9].  Nuckolls  [23]  used  a  value  of  1.0  .uega- 
bar  as  the  peak  shock  pressure  required  to  deposit  enough  heat  to  vaporize  the  tuff  In  his  cal¬ 
culation  of  the  phenomena  associated  with  the  RAliCIER  event. 

Bishop  handles  the  question  of  the  energy  available  to  the  shock  at  any  radius  by  deilnlrg 
an  "effective  shell  thickness."  The  c*scef«  is  based  on  the  fact  that  a  spherical  shell  carrying 
some  fraction  of  the  original  energy  of  !!«  explosion  eventually  breaks  arvay  from  the  hot  gas 
core,  and  once  this  has  happened,  no  energy  to  support  the  shock  can  be  received  from  the  core. 

In  the  strong  shock  region  with  which  the  discussion  is  presently  concerned,  the  shell  is  still  In 
complete  contact  with  the  core. 

Peet  [24,  p.  526]  points  out  that  Sharpe’s  calculations  of  the  seismic  signal  from  a  spherical 
cavity  show  that  the  shape  of  the  radial  pressure  profile  away  from  the  front  has  only  a  minor 
effect  on  the  results.  Both  theory  and  observations  indicate  that  the  total  pulse  duration  Is  short 
In  the  close-in  region  [23] .  Basic  communication  theory  asserts  that  If  100  cps  is  taken  as  the 
upper  limit  of  useful  energy  for  seismic  observations,  the  record  of  the  seismic  signal  from  a 
pulse  shorter  than  0.01  second  can  reveal  nothing  aboxd  the  wave  shape  or  spectral  content  of 
the  Input  [12,  25] .  Therefore,  no  recoverable  Inforitiation  Is  lost  If  the  shock  energy  In  the  source 
region  Is  considered  to  be  confined  to  a  thin  shell. 

The  effective  thickness  of  the  shell  Is  defined  as  having  a  volume  equal  to  tlie  total  energy 
In  the  shock  divided  by  the  energy  per  unit  volume  at  the  shock  front,  taken  to  be  the  kinetic 
energy  portion  of  Figure  1.  The  thickness,  y,  of  this  shell  In  the  strong  shock  region  Is  given 
by  Bishop  [13,  p.  48] ; 

In  which  the  subscript  f  refers  to  the  shock  front 
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and  R  is  the  radius  of  the  front.  The  radial  profile  of  density  p  and  particie  velocity  u  must  be 
known  to  apply  this  relationship. 

In  the  case  of  a  TNT  detonation  wav®,  y/R  ■  0.064,  which  means  the*  mopt  >  f  the  kinetic 
energy  Is  confined  to  a  sh»  ’!  bounded  by  the  front  and  an  inner  radius  only  6%  smaller  than  the 
ra'dius  of  the  front. 

Taylor's  theory  of  strong  shocks  f21]  provides  the  required  rad'.jl  variation  of  p  and  u. 
These  profiles  are  a  function  of  y,  the  adiabatic  exponent.  The  shell  thickness  is  t:.on  ^nly  a 
function  of  y,  which,  of  course,  is  a  property  of  the  medium. 

The  analysis  for  the  strong  shock  and  hydrodynamic  regions  will  be  summarized  on  the 
basis  of  Bishop's  procedure.  After  a  nuclear  explosion,  there  will  be  an  initial  zone  In  which 
the  parameter  x  >  100,  1‘aylor'r.  tVory  of  very  intense  shocks  gives  the  follov  ing  equations; 


-1  2 

Pj  =  2  (y  +  1)  p^c^  Xj 


(20) 


(21) 


where  P|  is  the  shock  front  pressure  for  x^  >  100  and  R^  is  the  radius  of  the  front  for  this 
pressure.  is  the  explosion  energy  release,  and  B(y)  Is  a  function  of  y  graphed  in  Figure  4. 
y  is  calculated  from  Equation  17. 


riOimE  4.  RELATIVE  SHOCK  PHESStmE  VS.  ADIABATIC  EXPONENT  (After  Bishop) 
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As  the  pressure  falls  below  the  "strong  shock"  level  (In  the  sense  that  Taylor's  equations 
no  longer  give  good  results),  but  is  still  hydrodynamic  (p^  >  Pj),  the  following  equations  may  be 
used  to  follow  the  shock  history: 


log 


fl  +5* 

*1  -1 
—  ♦  (2b  - 

X 

Dlog 

X 

1.-^ 

(22) 


(2S) 


where  p  =  overpressure  at  radius  R  In  the  hydrodynamic  range,  and  b  -  a  iieat  loss  parameter 
given  by 


b  -  1/3  (2  ♦  q^)  (24) 

where 

q^  -  0.4zj  (1  +  x^)  fl  -  1.82  (1  +  y)]’ *  (25) 

The  plastic  flow  limit,  x_,  Is  determined  from  experimental  Hugoniot  data  as  the  smallest 
m 

value  of  X  for  which  the  log  x  •  iog  p  graph  is  a  straight  line  (for  tuff  this  value  is  about  3.5; 

p 

see  Figure  3);  z  ■>  —  is  the  ratio  of  shock  front  radius  tc  effective  shell  thickness,  given  in 
Figure  5,  and  z^  is  the  value  at 


FIOimE  B.  SHELL  THICKNESS  VS.  ADIABATIC  EXPONENT  (Adrr  Bishop) 
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The  shell  thickness  ratio  In  the  hvdrodynamlc  region  is  given  uy 


Chaszeyka  f22]  describes  the  behavior  of  an  aggregate  of  silica  be^ds  and  air  intended  to 
simulate  some  actual  earth  materials  under  strong  shock  conditions.  When  a  shock  pr*^:;*is 
through  an  aggregate,  the  solid  particles  undergo  compression,  bending,  and  shear.  When  siz¬ 
able  voids  exist  in  the  aggregate,  the  compres.’ion  of  the  Included  air  causes  the  temperature 
of  the  mixture  to  increase  beyond  that  which  would  occur  if  a  uniform  sample  of  the  solid  were 
shocked.  The  resultliig  thermal  stresses  in  the  solid  add  to  the  destructive  effect  of  the  shock 
on  the  solid  component.  The  result  is  that  the  aggrogale  crushes  dow.i  to  its  solid  fraction  at  a 
pressure  that  is  low  compared  to  its  static  crushing  pressure. 

Even  though  the  shock  propagates  faster  in  a  gas  than  a  solid,  the  braking  effect  of  the  sur¬ 
rounding  solid  tends  to  hold  back  the  shock  through  the  air,  resulting  in  an  essentially  uniform 
shock  front. 

The  heat-transfer  process  from  the  compressed  air  to  the  solid  particles  is  complex,  in 
developing  the  Ranklne-Hugoniot  equation  of  state  for  the  aggregate,  Chaszeyka  assumes  that 
as  the  particles  are  crushed  down  (filling  the  voids),  they  break  up  Into  very  small  fragments 
which  are  thoroughly  intermixed  with  the  air  so  that  temperature  equilibrium  occurs  in  about 
one  microsecond.  His  equation  of  state  will  not  be  repeated  here,  but  the  resulting  cycle  for  a 
strong  shock  is  shown  In  Figure  6.  The  equatiov  of  state  Is  derived  by  shocking  the  air  and 
the  aggregate  to  the  peak  pressure,  and  then  letting  them  reach  an  equilibrium  temperature. 

Results  of  an  earlier  study  of  equations  of  state  and  blast-wave  properties  by  Chaszeyka 
and  Porzel  fl8]  are  shown  In  Flgiires  7  through  10.  The  materials  selected  were  two  types  of 
tuff,  sand,  clay,  sandstone,  granite,  and  ice.  The  overpressure-distance  curves  are  for  a  1-kt 
explosion. 

This  discussion  requires  an  additional  note  on  chemical  explosions  in  solids.  The  shock  is 
never  strong  (In  the  quantitative  sense  of  x  >  100  fl3,  p.  12B]),  and  it  is  known  f22,  p,  53-54] 
that  shock  waves  from  nuclear  explosions  do  not  scale  directly  to  those  produced  by  chemical 
oxplo.stons.  Bishop  explains  how  the  Initial  shell  thickness  ratio,  nt  the  boundary  Itotween 
the  detonation  wave  In  the  explosion  and  the  medium  can  be  calculated  In  terms  of  the  detonation 
velocity  in  the  explosive  and  the  pressure  at  the  boundary  between  explosive  and  medium.  This 
Input  pressure  may  not  bo  great  enough  to  make  the  process  hydrodynamic;  if  it  is  liigh  enough. 
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however,  then  It  is  a  straightforward  calculation  to  use  the  value  ot  to  get  the  other  Input 
quantities  needed  to  apply  Equations  22  through  25  for  the  hydrodynamic  region-  If  it  is  not,  the 
values  are  used  as  input  to  the  equations  In  the  cinishing  and  cracking  zone  to  he  developed  in  the 
next  section. 

Brode  [26]  provides  a  complete  example  of  the  calculation  o?  the  pressures  on  the  wall  of 
a  cavity  enclosing  a  chemical  explosion.  His  calculations  Indicate  that  the  cavity  pressure  is 
affected  by  the  loading  density  of  the  explosive,  the  equation  of  state  of  the  explosion  piodurts. 
the  presence  of  air  in  the  cavity,  and  the  details  of  the  detonation  of  the  explosive  material.  It 
is  of  greatest  importance  to  the  present  discussion  *hat  he  was  able  to  calculate  cavity  pressures 
which,  when  used  as  Inputs  to  the  calculation  of  elastic  waves  [27],  gave  results  agreeing  closely 
with  observed  values. 

Nicholls  [28]  has  also  done  recent  work  on  the  manner  in  which  the  properties  of  the  ex> 
plosive  affect  the  energy  appearing  in  the  form  of  seismic  waves.  He  concludes  that  the  maximum 
seismic  energy  is  generated  when  the  characteristic  impedance  of  the  explosive  (product  of  load¬ 
ing  density  and  detonation  velocity)  is  equal  to  the  Impedance  of  the  medium  (product  of  density 
and  eompressional  wave  velocity). 

All  investigators  agree  that  progress  in  the  analysis  of  strong  shocks  depends  on  additional 
experimental  data  on  the  behavior  of  earth  materials  at  the  pressures  involved.  As  we  have 
seen,  the  pressure-vs.-speclflc-volume  curve  over  the  entire  range  of  pressures  from  ambient 
to  the  highest  value  of  peak  shock  pressure  enters  Into  the  calculation  of  the  response  of  the 
medium. 

Cushing  and  Relly  [12]  point  out  the  fact  that  avaiaible  data  are  based  on;  (1)  static  tests, 

(2)  dynamic  tests  using  small  amplitude  elastic  waves,  and  (3)  dynamic  tests  using  shock  waves. 
Consequently,  there  is  a  lack  of  Information  about  the  behavior  of  soils  and  rocks  at  pressures 
In  the  neighborhood  of  the  crushing  strength.  The  result  is  that  the  response  in  the  terminal 
nonlinear  region  (discussed  In  the  following  section)  Is  the  greatest  source  of  uncertainty  in  the 
sequence  of  events  from  the  explosion  to  the  seismic  signal. 
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FIGURE  7.  PRESSURE- VOLUME  CURVE  FOR  TUFFS  (After  Chaereyka  and  Porr.el) 
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FIOtmE  S.  PEAK  OVERPnESSURE  VS.  DISTANCE  FOR  TUFFS  (1  klloton  yield) 
(After  Chnazeyka  and  Porael) 
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FIGURE  9.  PRESSURE-VOLUME  CURVES  FOR  VARIOUS  MATERIALS  (After  Chasr.eyka  and  Porzel) 
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PIOURK  10.  PEAK  OVERPRESSURE  VS.  DISTANCE  FOR  VARIOUS 
MATERIAIA  (1  klloton  yield)  (After  Chnsreykn  nnd  Porrel) 
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THE  TERMINAL  NONLINEAR  REGION;  THE  TRANSITION  FROM  A  SHOCK 
WAVE  TO  AN  ELASTIC  WAVE 

In  the  region  discussed  here,  the  material  no  longer  behaves  hydrodynatnicaily  because  the 
shear  strength  cannot  be  neglected,  but  the  stresses  still  exceed  th'»  elastic  limit.  The  processes 
taking  place  in  this  region  hold  the  key  to  many  of  the  most  ii.'.perfanf  questions  within  the  scope, 
of  this  report.  The  actual  transition  from  a  shock  wave  to  a  seismic  wave  occurs  at  a'-,  overpres¬ 
sure  less  than  the  crushing  strength  of  the  medium.  The  limits  of  this  region  of  transition  and 
the  characteristics  of  shock  propagation  within  It  are  controlled  by  the  crushing  stress,  the  yield 
stress,  and  the  elastic  limit.  These  properties  maybe  the  best  paramerers  to  use  in  describing 
quantitatively  the  effect  of  the  medium  on  the  seismic  output  of  an  explosion.  Furthermore,  the 
entire  pressure-volume  curve  for  the  p-  cssiire  range  corresponding  to  the  zon?  of  tra.isltion  will 
determine  the  waste  heat  deposited  in  It,  and  thereby  determines  the  energy  available  to  the  seis¬ 
mic  waves.  Many  investigators  retain  the  term  "pressure"  for  this  region,  using  it  synonymously 
with  "mean  stress,"  one  third  the  sum  of  the  principal  stresses. 

It  was  stated  above  that  there  Is  no  hope  of  recovering  Information  on  the  short-duration 
pulse  in  the  Initial  shock  from  seismic  waves  with  the  periods  usually  observed.  If  there  are 
processes  which  tend  to  lengthen  the  pulse,  however,  they  must  occur  in  this  zone  of  transition, 
and  it  may  be  possible  to  deduce  the  character  of  the  pulse  at  the  output  of  this  region  from  the 
recorded  elastic  waves  [12].  This  pulse  lengthening  is  not  an  apparent  increase  in  the  predomi¬ 
nant  period  because  of  selective  transmission  of  low  frequencies  by  the  earth. 

The  pressure  range  of  Interest  here  is  almost  unexplored  both  theoretically  and  experi¬ 
mentally,  except  for  work  on  metals  [13,  p.  69  J.  Ti.roretlcal  analyses  of  waves  in  this  region 
are  very  difficult  and  have  been  approached  through  the  application  of  theories  of  plastic  flow, 
fracture,  and  compaction. 

The  significance  of  this  region's  existence  has  long  been  recognized  by  exploration  seismol¬ 
ogists,  but  apparently  no  workers  in  this  field  have  seriously  attempted  to  determine  the  effects 
quantitatively,  Clewell  and  Simon  [29]  recognized  that  the  rise  time  of  the  input  seismic  signal 
is  very  much  lengthened  relative  to  the  detonation  pulse  from  a  high  explosive  because  of  inelastic 
processes  in  this  region.  The  result  is  a  limitation  on  the  power  of  reflected  waves  to  resolve 
thin  beds  in  the  sedimentary  section. 

Morris  [6]  also  discusses  the  types  of  effects  to  be  expected  in  this  region.  He  recognizes 
that  earth  materials  are  much  weaker  in  tension  than  in  compression,  and  points  out  that  the 
growth  of  cracks  will  accompany  tensile  stresses  and  that  energy  from  the  compressive  shock 
will  be  required  to  replace  the  strain  energy  released  through  cracking.  He  also  indicnte.s  the 
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change  in  pulse  shape  to  be  expected  because  lar^e  stresses  undergr  greater  attenuation  than 
small  ones  In  many  earth  media. 

This  last  point  is  made  clear  by  an  CKamlnation  of  the  dynamic  strcsi-e^rain  c..  **ve  for  silty 
clay,  Figure  11  [30],  Because  (iie  unloading  curve  is  almost  vertical,  considerable  compaction 
results  from  a  stress  cycle.  The  area  between  the  loading  {,r>d  unloading  curves  represents  en¬ 
ergy  lost  from  the  stress  wave.  Lampson  f30]  comments  thal  la-ge  stresses  not  only  undergo 
greater  attenuation,  but  also  propagate  with  a  lower  velocity,  so  that  the  waveform  is  spread 
out  in  time  and  diminished  in  amplitude.  The  seismic  wave  velocity  ordinarily  observed  is  de¬ 
termined  by  the  slope  of  the  curve  for  very  low  stresses. 

Peet,  in  discussing  seismic  waves  from  a  shock,  assumes  an  abrupt  transition  from  a  hydro- 
dynamic  region  to  elastic  behavior  f24|.  He  states  that  the  nonlinear  effects  may  t/C  confined  to 
a  small  region  between  the  two  zones,  or  the  region  may  scale  with  charge  weight  in  the  same 
way  the  shock  wave  region  does,  so  that  his  model  of  the  source  does  not  change.  The  boundary 
between  the  two  types  of  behavior  is  specified  in  terms  of  "a  certain  yield  stress,”  and  is  not 
further  defined.  Feet's  results  will  be  discussed  in  more  detail  in  the  section  on  elastic  wave 
propagation. 

Aokl  flO]  has  obtained  expressions  for  the  seismic  waves  generated  in  this  region  and  for  the 
wave  form  at  its  termination.  His  analysis  is  based  on  Tresca's  yield  condition  for  plastic  flow 
[31].  He  defines  a  plastic  wave  as  any  kind  of  wave  with  a  propagation  velocity  less  than  an 
elastic  wave.  But  since  these  Intermediate  stress  level  phenomena  may  propagate  in  the  form 


riounE  n.  expehimental  dynamic 

STHESS-STnAIN  CURVE  FOR  FREE  EARTH 
(Sn,TY  Ct^AY)  (After  Lampson) 
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of  cracks  produced  by  brittle  fracture,  Aoki  states  that  the  term  "fr-cture  wave”  may  be  more 

appropriate  than  "plastic  wave.”  He  considers  the  cases  of  waves  travellne  behind  the  elastic 

1/2  1/2 

wave  with  speeds  given  by  (k/p)  and  (p/p)  *■,  where  k  and  p  are  the  Incompressibility  and 
rigidity,  respectively.  In  either  ease,  an  elastic  precursor  runsaheadof  the  inain  shock,  and  as 
the  pressure  decreases,  the  two  merge  to  form  the  elastic  pulse. 

Cushing  and  Relly  (12]  have  established  an  elastic-plastic  static  mod.'.l,  using  a  procedure 
very  similar  to  Aoki's.  They  do,  however,  devote  more  attention  to  the  equation  of  state  i,t  the 
region.  They  have  also  started  work  on  a  d3mamic  model  In  which  the  previous  stress-strain 
history  of  the  medium  Is  taken  into  account,  an."  have  set  up  a  statistical  model  of  a  porous 
medium. 

Zvolinskii  (32]  assumes  a  iiiediv.m  in  which  there  is  an  abrupt  change  f-om  linear  elastic 
behavior  to  an  Incompressible  compacted  state.  Since  plastic  flow  In  the  compacted  state  is 
irreversible,  he  proposes  a  model  in  which  the  plastic  work  is  proportional  to  the  change  In  the 
greatest  shear,  with  the  proportionality  constant  depending  on  the  mean  stress.  Essentially,  he 
generalizes  the  work  of  Aokl  and  Cushing  and  Relly  by  using  a  more  general  yield  condition,  in 
which  the  difference  between  the  maximum  compresslonal  and  tensile  stresses  Is  a  function  of 
the  mean  stress  rather  than  a  constant.  He  is  one  of  several  authors  who  model  this  region  by 
allowing  the  material  to  crush  to  an  incompressible  state. 

Haskell  (33]  has  also  developed  a  static  theory  for  the  transition  zone,  using  an  incoherent 
granular  aggregate  rather  than  a  plastic  as  a  model  of  the  material.  Then  the  stress  In  the  zone 
of  failure  Is  determined  by  the  Coulomb- Mohr  .■'vlterlon  [31]  instead  ofaTresca  or  von  Mises 
yield  condition.  Haskell  compares  his  final  results  with  data  from  chemical  explosions  In  halite 
and  from  the  RAINIER  event.  He  finds  that  he  can  match  the  data  to  his  theory,  but  to  do  so  he 
must  use  values  for  the  internal  friction  parameter  that  are  much  lower  than  those  obtained  from 
compression  tests  on  unconsolidated  materials.  He  suggests  that  plastic  flow  and  fracturing  both 

4  * 

may  occur  In  the  zone,  giving  this  result. 

Fugelso  (11]  points  out  several  shortcomings  of  any  application  of  present  plastic  deforma¬ 
tion  theories  to  cases  In  which  the  loads  are  applied  very  rapidly,  or  for  a  very  short  time.  The 
Inherent  difficulties  In  developing  a  theory  of  dynamic  plastic  deformation  result  from  the  nonlin¬ 
earity  of  the  equations  and  the  fact  that  the  stress-strain  curve  Itself  Is  a  function  of  the  loading 
rate.  Fugelso  has  developed  a  theory  of  plastic  deformation  based  on  the  theory  of  dislocations 
within  a  crystalline  body,  and  applied  It  to  the  problem  of  Impact  loading  of  a  rod,  combining  the 
effects  of  elastic  distortion  and  the  movement  of  dislocations. 
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Fugelso's  conclusions  Include  the  following: 

1.  The  deformation  resulting  from  the  ar-iva!  of  the  stress  wave  3  -urri  in  two  steps:  an 
instantaneous  elastic  deformation  and  a  time -dependent  plastic  defu  'nnMnr  . 

2.  The  elastic  wave  propagates  a«  a  distinct  wave,  but  ther^  no  distinct  plastic  wave. 

3.  The  stress-strain  curve  for  any  material  is  not  uiliq'in,  but  is  dependent  on  the  strain 
rate  and  the  duration  and  magnitude  of  the  stress. 

In  view  of  the  work  done  also  by  Aoki  and  Zvolinskli  (above),  the  conclusion  reached  here  is 
important:  that  th**  solution  gives  no  forward  uropagatlng  wave  that  moves  with  a  velocity  lower 
than  that  of  the  elastic  compressional  wave.  The  solution  is  valid,  of  course,  only  for  the  one 
type  of  nonelastic  deformation  mechanism  assumed.  Pugelso  states  that  his  mode’,  is  intended 
for  crystalline  bodies,  and  soils  are  specifically  mentioned  among  n.'dia  for  which  more  com¬ 
plicated  mechanisms  may  be  involved. 

An  examination  of  the  pressure -volume  diagram  for  a  soil-like  material  (see  Figure  12)  in 
terms  of  Equation  4  indicates  part  of  the  reason  for  the  complicated  nature  of  the  transitional 
processes  f  12,  17,  19],  The  velocity  of  shock  propagation  is  proportional  to  square  root  of 
(p  -  pQ)/(Vg  -  v).  The  linear  portion  of  the  pressure-volume  curve,  with  p^  as  the  elastic  limit, 
fixes  a  slope  which  divides  the  curve  into  two  parts.  For  a  pressure  Pj,  (p^  -  Pq)/(Vq  -  v^)  > 

(p^  -  Pq)/(Vq  -  Vj),  the  shock  is  ’’superselsmlc,"  and  the  conditions  for  a  single  stable  shock  are 
satisfied  fl7].  However,  for  a  pressure  Pg  below  the  intersection  of  the  elastic  portion  (extended) 
with  the  pressure -volume  curve,  the  InequalPy  Is  reversed  and  the  shock  breaks  up  into  a  two- 
wave  structure.  Figure  13,  taken  from  CushiiA  .and  Relly  [12,  p.l6],  shows  the  radial  profile  of 
the  resulting  pressure  pulse.  The  arrival  of  the  front,  traveling  at  the  seismic  compressional- 
wave  velocity,  is  accompanied  by  a  steep  rise  in  stress  to  the  elastic  limit.  The  pressure  then 
increases  gradually  to  an  Intermediate  value,  Pj,  determined  by  the  tangent  to  the  equatlon-of- 
state  curve  through  Pg,  Vg.  Relative  to  the  p^,  v^  state,  conditions  for  a  stable  shock  exist,  so 
that  the  steep  rise  to  the  peak  pressure  Pg  occurs. 

Some  important  consequences  of  this  behavior  are  (1)  a  zone  for  which  the  peak  stress  is 
well  above  the  crushing  strength,  but  for  which  the  peak  stress  moves  subseismically,  so  that 
elastic  precursors  result;  and  (2)  a  mechanism  for  lengthening  the  pulse  in  time.  A  computed 
pressure-time  history  for  a  point  in  the  transition  zone  of  the  RAINIER  event  is  seen  in 
Figure  14  [23]. 

Chaszeyka  [22]  conducted  a  systematic  Investigation  of  the  first  stages  of  the  transition 
process  both  In  hard,  low-poroslty  material,  and  in  dry,  unconsolidated,  noncoheslve  soils,  such 
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FIGURE  14.  COMPUTED  STREM  WAVE  IN  NONLINEAR  REGION,  RAINER 


EVENT  (After  Nuckolls) 

as  dry  sand.  His  study  partially  supports  the  assumption  made  by  Peet:  nonelastic  processes 
do  not  greatly  affect  the  generation  of  seismic  waves  in  the  case  of  a  porous,  sandy>type  soil. 
For  a  granite-like  rock,  however,  a  seismic  wave  of  considerable  amplitude  is  generated  when 
the  pressure  in  the  shock  front  is  at  or  above  the  Trushing  strength.  The  zone  of  crushing  and 
cfcmpaction  Is  examined  in  this  work. 

For  a  solid  which  is  malleable  and  nearly  isotropic  and  has  a  negligible  air  content,  the 
shock  weakens  through  gradual  decrease  of  the  deformation  until  it  becomes  an  elastic  wave. 
The  transition  is  not  sharp.  If  the  solid  is  porous  or  if  it  is  equivalent  to  an  aggregate  of  solid 
particles  and  air,  the  propagation  changes  considerably  during  the  transition. 

A  porous,  aggregate-type  medium  will  have  a  lower  crushing  limit,  p^j,  at  which  the  solid 
grains  must  begin  to  fracture,  and  an  upper  crushing  limit,  p^^,  at  which  the  interstices  are 
filled  with  crushed  material  (see  Figure  15).  If  hydrodynamic  behavior  is  defined  as  occurring 
above  the  crushing  strength,  the  transition  begins  near  p  .  Between  p  and  p  ,  is  a  region  of 
•  invariant  reduced  volume  (Chaszeyka's  terminology),  or  compaction.  The  pressure  in  the  shock 

drops  from  p^^  because  of  thermal  losses  and  work  on  the  medium.  At  a  pressure  p^  between 
the  two  crushing  limits,  the  volume  Is  reduced  to  a  fixed  value  v^.  This  behavior  is  indicated 
by  the  vertical  dashed  lines  on  Figure  15.  The  transition  from  compacting  to  plastic  liehavlor 
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occurs  near  p^j;  the  final  transition  to  a  seismic  pulse  occurs  at  an  overpressure  below  the 
limit  for  plastic  flow. 

The  solid  particles  are  compressed  durln({  the  compacting  process,  bul  ‘lie  r lunge  in  volume 

from  this  cause  is  negligible  compar®/!  to  that  arising  from  fracturing  of  the  grains.  The  air 

shock  in  this  region  Is  low.  The  energy  loss  that  attenuiteh  the  shock  from  p  to  p  ,  is  nearly 

cu  *^cl  ^ 

all  due  to  work  done  in  crushing  the  solid  portion  of  the  aggregk.'e. 

The  behavior  of  a  hard  low-poroslty  rock  Is  indicated  by  the  step  on  Figure  15.  Hire  a  de¬ 
finite  value  of  the  crushing  stress  exists  at  which  the  material  crushes  to  hydrodynamic  behavior; 
below  this  value  it  behaves  nearly  elastically. 

Ctiaszeyka  gives  a  form  for  the  variation  in  density,  particle  velocity,  and  pressure  through 
the  transition  aone.  Since  the  density  to  which  the  material  is  compacted  decreases  as  the  peak 
pressure  drops  with  increasing  radiki  ^distance,  the  density  will  be  low^r  at  the  shock  front  than 
behind  It.  A  density-distance  relation  of  the  form 


P  =  Pg(R^r)q  (27) 

is  assumed,  in  which  r  is  the  distance  to  a  point  behind  the  front  (r  <  R),  and  q  =  - 1  P-. 

d  r 

By  conservation  of  mass  (and  using  the  fact  that  p  is  constant  behind  the  shock  in  the  com¬ 
pacted  zone),  the  following  expression  for  q  is  obtained; 


q  e  3 


p,  1  - 

pQ  1  -  (R/R)^ 


(28) 


in  which  R^  is  the  radius  corresponding  to  p^^,  and  is  the  smallest  radius  at  which  compaction 
to  constant  density  occurs. 

For  a  particle  velocity  u^  at  the  front,  the  velocity  at  a  point  in  the  zone  is 


u  ■ 


UsP.R- 


(29) 
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The  pressure  distribution  is  found  from  conservation  of  momentur.<.  In  terms  of  quantities 
already  defined, 

^  Pi/f’o  1  l^d  log  (p  -  d  log  R  j' 

This  equation  shows  that  for  an  aggregate  the  pressure  decreases  toward  the  origin  from  the 
shock  front. 

The  decrease  of  peak  pressure  with  radius  of  the  shock  is  found  by  r  technique  ol  blast  wave 

integration  developed  by  Chaszeyka  and  Porzel  fl8].  The  computed  result  for  an  aggregate  of 

“6 

69%  silica  and  31%  air  Is  shown  in  Figure  16;  p  falls  off  at  about  R  .  Comparison  with 
Figures  2  and  3  shows  that  the  pressure  range  covers  the  region  of  the  crushing  zone  designated 
by  Bishop.  The  radius  at  which  the  upper  crushing  limit  occurs  can  be  found  for  any  material 
from  a  curve  like  the  one  shown  In  Figure  16  if  adequate  experimental  information  is  available. 

For  a  porous  medium  in  the  zone  R  >  R^,  the  value  of  q  in  Equation  27  varies  from  about 
1.6  (when  R  is  10%  greater  than  R^)  to  zero  (as  R  approaches  infinity).  In  this  case,  the  small 
part  of  the  original  energy  from  the  explosion  still  remaining  (less  than  3%)  is  dissipated  in 
crushing  the  solid  particles;  the  shock  weakens  as  a  result  of  both  this  crushing  work  and  spher¬ 
ical  spreading. 

For  a  hard  rock,  approaches  rapidly  as  the  stress  drops  below  the  sharply  defined 
crushing  limit.  The  exponent  q  has  a  low  value  and  approaches  zero  quickly,  being  about  0.02 
when  R  is  10%  greater  than  R^.  The  result  is  a  relatively  small  energy  decrease,  and  a  strong 
seismic  pulse  will  be  generated  at  a  pressure  only  slightly  below  the  crushing  limit. 

The  transition  discussed  here  Is  only  the  first  in  the  sequence  leading  to  the  seismic  pulse 
[22,  p.  144].  For  later  transitions,  associated  with  lower  overpressures,  we  return  to  Bishop's 
treatment.  He  divided  the  equation  of  state  into  five  regions,  with  several  subregions  (see 
Figures  2  and  3).  Starting  from  the  high  pressure  end.  Region  V  is  the  previously  discussed 
hydrodynamic  region  for  which  x  >  1,  and  y  is  nearly  constant.  A  subrange  in  which  large  plastic 
deformations  occur  is  seen  Just  above  the  hydrodynamic  termination  pressure,  Ph- 

The  pressure  Pjj  is  approximated  by  Pj,  corresponding  to  x  0.5,  and  marks  the  upper  limit 
of  Region  IV,  the  zone  of  crushing.  The  crushing  pressure  p^  determines  the  lower  boundary  of 


P-Po“^Ps-Po> 
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crushing.  The  crack  zone,  Region  Til,  is  subdivided  into  two  parts:  IIIp,  a  region  of  hoop-stress 
tensile  cracks,  mb,  in  which  shear  stress  cracks  cccur.  The  lower  limit  of  the  crack  zone  is 
designated  p^,  which  denotes  the  highest  pressure  In  tue  stress  wave  for  whlcn  processes  are 
reversible.  The  boundary  between  the  two  types  of  crack  occurs  at  the  maxlmutiu  In  the  x  -  p 
curve,  at  a  pressure  near  twice  the  ambient  pressure.  Regions  HI  cv't  IV  represent  irreversi¬ 
ble,  nonhydrodynamic  behavior. 


Region  n  is  that  in  which  processes  are  reversible,  but  strains  are  finite;  Heglon  i  )s  the 
zone  of  elastic  waves.  Region  n  (called  the  "semiacoustic  zone"  by  Bishop)  begins  at  an  over¬ 
pressure  of  about  0.02  p^. 

The  low-pressure  end  of  the  x  -  p  diagram  will  be  considered  In  detail  first.  For  strictly 
infinitesimal  strains,  the  adtabaUc  exponent  is  =  ^o®oV*’o  derivation  of  Equation.  13).  For 
condensed  media  (liquids  or  solids),  y^  is  large  compared  to  one  for  :  ear  one  bar.  For  small 
but  finite  strains,  the  transition  across  tht?  shock  will  be  nearly  reversible.  As  Porzel  [19,  p.  41] 
states,  all  materials  are  excellent  transmitters  of  shock  energy  at  low  pressures.  Bishop  postu¬ 
lates  a  form  for  y  for  the  finite-strain,  nearly  reversible  case  that  Is  analogous  to  that  for  y^ 
above: 


y 

'rev 


(31) 


in  which  y^^^  is  the  semlacoustic  exponent,  and  the  compresslonal  wave  velocity  c^  has  been 
replaced  by  the  shock  velocity  U.  y^.^^  is  also  a  large  number,  so  that  the  value  of  x  for  reversi¬ 
ble  processes  is  approximated  by 

x^  =  1/2  'x/p^)  (32) 

where  p  is  small  enough  for  the  assumption  of  reversibility  to  hold.  It  follows  that  x^  must  be 
small  compared  to  one.  In  order  to  carry  on  without  experimental  data.  Bishop  assumes  O.S  as 
the  maximum  value  of  x  for  the  semiacoustic  region.  From  Equation  32,  this  requires  that  p^, 
the  reversible  limit  in  terms  of  pressure, be  less  than  or  equal  to  p^,  the  ambient  pressure. 

This  choice  of  the  reversible  limit  is  based  on  the  fact  that  most  earth  materials  have  very 
low  tensile  strengths.  They  will  fail  in  tension  as  soon  as  the  tensile  hoop  stresses  associated 
with  the  radial  compressive  stress  overcome  the  llthostatlc  compression.  This  will  occur  at  an 
overpressure  close  to  the  ambient  pressure. 

By  analogy  with  the  elastic  case,  dynamic  bulk  and  shear  moduli  can  be  defined  such  that 
K  *  K  (1  I  x)  «  k  +  4/3  ps.  In  the  semlacoustic  region,  it  is  not  unreasonable  to  take  k  /^a 

9  Q  9  S 

■  k/p,  and  Bishop  postulates  that  this  relation  may  be  approximately  true  up  to  p  =  2p^,  the  maxi¬ 
mum  point  shown  on  Figures  2  and  3. 
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The  behavior  In  crushing  has  been  discussed  In  the  review  of  Chaszeyka's  work;  Bishop's 
treatment,  however,  will  be  summarized  here  for  the  sake  of  continuity,  tn  the  crushing  region, 
the  dynamic  bulk  modulus  k^  is  almost  constant,  v  hlle  the  dynamic  shea:*  modulus  ps  Is  very 
small.  The  x  -  p  curve  Is  pr-:ir*iy  fiat,  with  a  constant  value  of  x  given  by 

x^  =  -1  +  2k/K^  =  1  -  (33) 

(see  the  discussion  following  Equation  16).  Bishop  concludes  that  no  equation  of  state  can  be 
established  for  the  crushing  zone  except  by  measurements  on  specific  materials. 

The  onset  of  Region  III  (coming  up  from  ihe  low  pressure  end)  occurs  when  the  medium  is 
first  broken  up  into  large  fragments  by  radial  cracks  from  the  hoop  stresses.  This  occurs 
when  the  tensile  stress  exceeds  the  remiltant  of  the  ambient  stress  and  the  tensile  strc.ig'h  of 
the  medium.  It  is  not  possible  to  assign  a  precise  value  to  at  which  cracking  begins  because 
of  the  imperfect  nature  of  earth  materials.  The  value  p^  seems  io  be  a  reasonable  estimate. 

As  the  stress  increases  above  p^,  the  fragments  are  broken  into  smaller  and  smaller  pieces 
until  the  medium  is  reduced  to  an  almost  homogeneous  aggregate  of  particles  of  relatively  high 
strength;  crushing  then  begins  at  pressure  p^.  Radial  cracking  predominates  up  to  a  pressure 
of  about  ?p^;  above  this  pressure,  circumferential  cracks  produced  by  shear  stresses  combine 
with  the  radial  cracks  to  make  the  smaller  pieces.  In  the  zone  of  shear  stress  cracks  (see 
Figures  2  and  3)  the  x  -  p  curve  has  a  negative  slope,  the  shock  velocity  decreases  with  increas¬ 
ing  pressure,  and  the  previously  discussed  elastic  precursors  result. 

Bishop  has  developed  equations,  on  the  basic  of  numerous  assumptions  and  approximations 
necessitated  by  the  lack  of  experimental  data,  •■'hich  make  It  possible  to  follow  the  history  of  the 
shock  through  the  transition  zone.  He  discusses  these  equations  thoroughly;  no  attempt  to  Justify 
them  fully  will  be  made  here. 


The  output  of  the  hydrodynamic  zone  is  the  input  to  the  crushing  zone.  Therefore,  the  cal¬ 
culations  start  by  finding  Rj  and  Pj,  the  approximate  radius  and  pressure  of  the  shock  at  hydro 
dynamic  termination,  from  Equations  22  through  25.  The  effective  shell  thickness  ratio  Zj  Is 
found  from  Equation  26.  Then,  for  p^  <  p  <  Pgt 


log^  =  (l+Xc)'*log.^ 


(34) 


»  I  (1 


+  1.5x^)  log- 


*JVc 


(35) 
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where  is  given  by  Equation  33. 

Some  assumption  about  losses  from  the  shock  as  It  passes  over  the  mr'iium  is  essential  to 
the  development  of  these  equations.  Th<=  blaR*.  wavs  integration  procedure  ilcvo’oppri  >  y 
Chaszeyka  and  Porzel  fl8]  Is  «»  "lethod  of  keeping  track  of  the  waste  heat  deposited  during  prop¬ 
agation  of  the  shock.  Bishop  treats  the  problem  In  the  followit manner. 

The  heat  loss  per  unit  mass  from  the  shell  Q  is  found  In  terms  Of  the  waste  heat  Q  (defined 
above)  by  =  G^G^Q,  In  which  G^  is  a  shock-strength  factor  which  is  zero  for  very  strr-.g 
shocks,  and  G^  is  a  factor  which  depends  on  the  medium.  By  assuming  that  G^  =  (1  +  x)'^,  it 
Is  possible  to  write  G^  -  (1  +  x)  Q^/Q-  A  value  of  G  =2.0  was  found  to  agree  well  with  experi¬ 
mental  data  in  the  crushing  zone  for  halite,  tuff,  and  alluvium,  the  only  earth  materials  investi¬ 
gated.  Equations  34  and  36  are  based  on  this  value. 

In  the  cracking  zone,  with  the  overpressure  between  the  crushing  and  reversible  pressure, 
good  results  were  obtained  by  the  use  of  arj  Inverse  square  law  for  ot  erpressure  decay; 

log  (p^/p)  =  2  log  (R/R^)  (36) 

where  R  Is  the  shock  radius  at  which  p  occurs.  '  - 
c  c 

In  the  semi-acoustic  zone,  p  <  p^. 


The  results  of  applying  the  theory  to  available  data  from  explosions  in  halite  (COWBOY  and 
GNOME  events)  and  volcanic  tuff  (HOBO  and  RAINIER  projects)  are  shown  in  Figures  17  and  18. 
Equatlon-of-state  data  were  taken  from  Lombard  [34],  and  ambient  stress  values  are  taken  from 
Adams  and  Swift  [35].  The  COWBOY  data  are  taken  from  Murphey  [36];  the  GNOME  values  are 
from  preliminary  results  obtained  by  the  Sandia  Corporation. 

On  these  figures,  R,  Is  the  radius  In  feet  of  a  sphere  of  pelletol  TNT.  With  the  pelletol 
•  S  -1/3 

density  taken  as  1.0  gm/cm  ,  R/Rj  =  6.38  R  W  '  ,  where  W  is  the  charge  weight  in  pounds. 

The  value  of  R^  for  the  GNOME  npelear  event  was  taken  ao  28.5  ft,  the  radius  of  a  pelletol  sphere 
with  the  same  energy  release,  S  kllotons.  The  cube-rbdt  scaling  brings  the  nuclear  and  small 
chemical  explosion  data  into  good  agreement.  » 

Using  the  properties  of  pelletol  and  halite.  Bishop  sihows  that  the  input  pressure  of  100 
kllobars  corresponds  to  the  upper  limit  of  the  zone  of  crushing.  There  is  consequently  no  hydro- 
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dynamic  region,  and  all  pressures  from  the  COWBOY  series  arc  In  *he  cracking  region.  The 
plotted  pressures  In  the  figures  were  not  measured,  hut  were  calculated  from  observed  particle 
velocities,  using  equations  developed  In  the  theory. 

For  tuff,  the  Input  pressure  for  pelletol  Is  75  kllobars.  Values  derived  from  Equations  34 
and  35  are  shown  on  Figure  18.  The  first  quantitative  results  On  the  effect  of  source  depth  en¬ 
countered  In  this  review  can  be  seen  In  the  HOBO  data.  The  rC'  ^-slble  limit  p  Is  seen  to  in- 

r 

crease  significantly  as  the  depth  of  the  shot  (and  therefore  the  ambient  stress)  Inc-ea^es.  The 
maximum  value  this  limit  can  attain  has  been  set  as  the  resultant  of  the  ambient  stress  and  the 
tensile  strength  of  the  medium.  The  actual  value  may  be  smaller  than  this  because  of  pre¬ 
existing  cracks  or  other  departures  from  uniformity  of  the  medium.  For  the  tuff,  the  reversible 
limit  Is  near  the  maximum,  hut  It  appears  to  be  much  less  than  the  maximum  of  59  birs  fc;-  the 
COWBOY  data  on  Figure  17. 

The  two  methods  of  analysis,  Cfiaszeyka's  and  Bishop's,  can  be  compared  by  the  results 
presented  on  Figures  8  and  18  for  explosions  in  tuff.  The  abscissas  of  Figure  8  are  to  be  divided 
by  6.0  (the  value  of  In  meters)  for  one  kiloton,  and  the  logarithm  taken.  Thus  six  meters  on 
Figure  8  corresponds  to  0.0  on  Figure  18,  30  meters  to  0.7,  and  60  meters  to  1.0.  The  region  of 
overlap  Is  therefore  limited,  but  within  this  region  the  agreement  Is  very  good.  The  tendency  for 
the  pressure-distance  curve  to  be  slightly  concave  to  the  distance  axis  shown  on  Figure  18  is  not 
seen  on  Figure  8.  It  Is  concluded  that  both  methods  of  analysis  are  sound,  and  the  principal  need 
is  for  more  experimental  data. 

Although  the  nonlinear  processes  have  been  analyzed  for  the  manner  In  which  they  absorb 
energy  from  the  shock  pulse,  no  analysts  has  be  >>;  made  of  the  seismic  signal  proceeding  from 
these  disturbances  directly.  Only  the  radial  compresslonal  pulse  Is  followed.  A  growing  crack 
Is  a  source  of  seismic  waves,  and  experimental  evidence  has  been  accumulated  that  these  waves 
can  be  observed  [37].  It  has  been  suggested  also  that  the  boundary  between  the  expanding  explo¬ 
sion  products  and  the  surrounding  medium  may  be  unstable,  developing  wave-llke  perturbations 
SO  that  the  output  from  the  source  may  not  be  spherically  symmetrical  [38] . 

Propagation  In  the  elastic  region  will  be  reviewed  next  on  the  assumption  that  the  source  Is 
spherically  symmetrical;  the  evidence  of  asymmetry  will  be  discussed  subsequently, 

6 

THe  IlASTIC-WAVe  RtOION 

The  entire  earth,  outside  of  the  hydrodynamic  and  nonlinear  regions  discussed  In  the  pre¬ 
vious  sections,  conslllutos  the  region  of  elastic  wave  propagation.  This  report  does  not  consider 
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the  details  of  the  effects  of  departures  from  Idea!  elasticity  in  this  region  on  the  shape  and  am¬ 
plitude  of  the  seismic  signal.  Also  omitted  are  the  derails  of  the  structure  of  ihe  crust  and  upper 
mantle,  which  must  be  taken  into  account  in  deter nilning  the  actual  amplitude  of  the  various 
seismic  waves  as  observed  at  dist^  it  points  fl].  Attention  is  directed  to  the  signal  at  short 
ranges,  and  at  long  ranges  if  the  medium  is  homogeneous  and  in|!i.:te  in  extent. 

Theoretical  and  experimental  studies  of  explosion-generated  elastic  wo.vos  have  been  car¬ 
ried  out  for  many  years,  starting  with  the  early  work  In  exploration  seismology,  sofuewhat 
later  extending  to  quarry  and  mine-blast  observations,  and  most  recently  to  nuclear  explosions. 
Some  of  the  experimental  work  has  been  designed  'o  test  theoretical  results,  but  most  of  it  has 
been  directed  to  empirical  correlations  of  maximum  seismic  wave  amplitude  with  various 
parameters  characterizing  the  explosion  and  observation  conditions.  Results  o!  the  search  tor 
empirical  scaling  relations  are  covered  io  a  later  section  of  this  report.  The  early  theoretical 
papers  made  no  attempt  to  determine  numerical  values  of  amplitude  in  terms  of  charge  weight. 

In  a  classic  paper,  Sharpe  [3]  attempted  to  explain  certain  facts  known  from  experience  to 
exploration  seismologists.  Referring  specifically  to  seismic  reflections  from  horizons  in  the 
sedimentary  section,  he  summarized  these  observations  as  follows: 

1.  A  given  amount  of  explosive  detonated  in  a  clay  or  water -saturated  sand  formation  re¬ 
sults  in  a  greater  amplitude  of  reflected  motion  than  an  equal  charge  detonated  in  a  dense,  rigid 
formation  such  as  limestone. 

2.  If  a  hole  is  sprung  by  an  initial  large  charge  in  order  to  form  a  sizable  cavity,  later 
small  charges  will  result  in  a  larger  amplitude  of  ,'iiflccted  motion  than  would  be  produced  In 
the  absence  of  springing. 

3.  The  frequency  spectrum  of  reflected  motion  is  a  function  of  the  formation  in  which  the 
charge  is  fired:  shots  fired  In  the  low  velocity  zone  result  in  very  low  frequency  motion  com¬ 
pared  to  shots  fired  below  the  low  velocity  zone;  shots  fired  in  a  rigid  material  (such  as  a 
limestone)  result  in  a  much  higher  frequency  motion  than  shots  made  in,  for  example,  a  shale; 

In  a  general  way,  the  high  frequency  content  of  reflected  motion  Increases  with  an  increase  in 
shooting  depths. 

4.  The  frequency  spectrum  of  reflected  motion  is  a  function  of  charge  size;  a  large  charge 
has  a  tendency  to  increase  the  proportion  of  low  frequencies  in  the  reflected  motion. 

5.  The  amplitude  of  reflected  motion  produced  by  a  given  quantity  of  a  high  explosive  Is 
much  greater  than  that  produced  by  a  quantity  of  low-speed  explosive,  even  when  the  maximum 
pressure  is  the  same. 
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These  statements  touch  upon  every  point  of  concern  within  this  pres..nt  report;  i,e.,  the  effect 
on  the  amplitude  and  spectral  content  of  the  seismic  clgital  of  the  medium  S'^ound  the  shot,  the 
depth  of  burial,  the  yield,  the  coupling,  and  tlie  properties  of  the  explosive  material. 

6.1.  THE  WAVEFORM  OF  EXPLOSION-GENERATED  P  wave? 

Sharpe  took  as  his  source  model  a  spherical  cartty  ariurd  the  shot  having  a  radius,  a,  iserge 
enough  so  that  only  elastic  processes  took  place  at  or  outside  the  surface.  This  cavity  is  desig¬ 
nated  the  "equivalent  cavity"  or  the  "equivalent  radiator"  by  various  authors.  An  exponentially 
decaying  pressure  pulse,  with  the  pressure  equal  at  all  points  of  the  spherical  surface,  was 
postulated  as  the  loading  equivalent  to  the  explosion;  p{a,  t)  = 

No  attempt  was  made  to  relate  a,  p^,  or  o  to  the  size  of  the  charge  and  the  pror-erlies  of 
the  medium. 

Although  Sharpe  obtained  his  results  for  the  special  case  in  which  Poisson's  ratio  a  >  0.2R. 
or  X  e  4,  it  is  useful  to  generalize  the  result  for  an  arbitrary  elastic  medium.  This  has  beei; 
done  by  Blake  [5],  but  the  solution  will  be  given  in  a  slightly  different  notation.  The  ratio  of  shear 
modulus  to  Lamp's  constant,  m/X  »  m,  will  be  used  In  place  of  the  parameter  K  introduced  by 

Blake.  Thus,  a  ■  JTfVjnl*  m  ■  1  for  o  »  0.25. 

The  displacement  potential  for  the  exponential  pressure  on  the  cavity  wall  is  given  by 


in  which  p  •>  density 


r  ■  distance  from  center  of  cavity 


2cVtn  +  m^ 
o  "  TTi' 


compressional  wave  velocity  in  the  medium 


(38) 
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The  displacement,  which  is  purely  radial  under  the  assumptions,  .s  the  derivative  of  the 
potential  in  Equation  38.  For  the  important  special  case  of  a  step  in  pressure,  0=0,  and  the 
displacement*  is 


The  first  term  represents  the  permanent  displacement  resulting  from  the  fact  that  the 
pressure  p^  continues  to  act  on  the  cavity  forever.  At  large  ranges  (r  >>  a),  the  last  term  in 
Equation  39  predominates. 


The  displacement  is  thus  a  damped  sinusoid  with  angular  frequency  o)^,  which  varies  directly  as 
the  compresslonal  wave  velocity  and  inversely  as  cavity  radius,  and  becomes  low  for  small 
values  of  m,  characteristic  of  soils. 

For  small  values  of  m,  the  amplitude  becomec  large  (and  4  becomes  small),  assuming  that 
the  peak  pressure  and  radius  at  the  inner  bounda^  %  oi  the  elastic  region  are  unchanged. 

The  damping  factor  (fraction  of  critical  damping)  C  is  &  function  of  the  properties  of  the 
medium; 


The  motion  in  rock  (for  which  m  Is'neair  1.0)  will  be  more  highly  damped,  or  less  oscillatory, 
than  In  soil,  Aokl  [10]  presents  waveforms  for  particle  velocity  which  Illustrate  this  for  m  =  1.0 
and  1/75  (see  Figures  19  and  2.0).  The  effect  of  including  elastic  precursors  from  the  nonlinear 
region,  as  Aokl  does,  is  to  increase  the  rise  time  of  the  signal.  The  following  notation  applies 
to  these  two  figures; 

*The  syinlx)!  "u"  Is  widely  used  for  particle  velocity  in  hydrodynamics  and  particle  displace¬ 
ment  in  elasticity  theory.  Rather  than  Introduce  a  new  symbol,  the  same  symbol  is  used  for 
both  in  different  parts  of  this  report.  From  this  point,  "u"  is  displacement,  "v"  Is  velocity. 
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V  =  particle  velocity 

b  =  radius  of  beginning  elastic  zone  a) 
k  =  constant  from  von  Mlses*  yield  condition:  o 


r-b 


rr 


T  =  — — ^ - : — .  t  is  time  at  which  shock  reaches  the  radial  distance  b 

b  b  '  o 

ac  =  propagation  velocity  of  "plastic"  or  "fracture"  wavc;  tiiken  as-v  —  ,  where  E  Is  the  bulk 
modulus  for  "plastic"  waves  and  the  rigidity  modulus  fcr  "fracture  "  waves 
c  *  compresslonal  wave  velocity 

A.  /V 

s  *  decay  exponent  in  pressure  pulse  in  notation  of  this  report) 

In  Figure  19,  w’lerc  A  =  p,  m  =  1,  the  two  velocities  of  inelastic  wave  propagation  are 
and  for  the  "plastic"  and  "fracture"  waves,  respectively.  In  Figure  19,  s  -  0,  cor¬ 
responding  to  a  step  function  in  pressure. 


The  results  for  a  clay-like  material  (>  *  75m;  see  Flgura  20)  are  baaed  on  a  velocity  of  the 
elastic-plastic  boundary  of  0.9  c  and  four  rates  of  pressure  decay.  The  initial  spike  (corre¬ 
sponding  to  rapid  pressure  decay,  or  large  values  of  s)  will  be  highly  attenuated  at  a  distance 
because  of  the  transmission  characteristics  of  the  earth,  discussed  below. 


The  response  to  arbitrary  forcing  functions  can  be  obtained  by  superposition  of  solutions 
of  the  forms  given  above.  Duvall  [4],  working  with  radial  strain  data,  was  able  to  arrive  at 
waveforms  closely  resembling  his  field  records  by  adding  the  responses  to  two  exponential 
pressure  functions  with  different  rates  of  decay.  He  found  that  at  short  ranges  the  non-oscllla- 


^ 
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tory  part  of  the  solution  predominates,  and  the  time  of  the  first  zero  crossing  Is  not  accurately 
predicted  by  taking  a  half  period  corresponding  to  He  also  points  out  that  the  range  at  which 
elastic  behavior  begins  cannot  be  determined  bv  locating  the  limit  of  non-oselllalurv  wave  shapes. 

Van^  [39,  40,  41]  has  c^^tntned  »t»e  effect  of  the  higher  power  terms  In  1/r  In  solutions  like 
Equation  39,  These  terms  produce  a  region  two  or  three  tin.cs  the  radius  a  In  which  the  peak 
amplitude  of  the  spherical  compresslonal  waves  decays  faster  Oan  r"^'  Since  the  material  is 
assumed  to  respond  In  a  perfect  clastic  manner,  this  requires  an  abnormal  attenuation  due  to 
geometric  spreading.  Van^  concludes  that  the  peak  amplitude  propagates  faster  than  the  wave 
front  In  this  zone,  the  more  rapid  propagation  yielding  the  faster  attenuation.  This  conclusion 
implies  that  the  beginning  of  the  elastic  zone  cannot  be  determined  experimentally  by  finding 
the  shortest  range  at  »  hl(  h  the  peak  displacement  decays  as  the  Inverse  first  powo.  of  the  dis¬ 
tance. 

An  excellent  upporbinlly  to  test  thl^  entire  approach  to  the  theory  of  explosion-generated 
seismic  waves  was  offered  by  the  Project  COWBOY  experimental  series.  The  objective  of  the 
program  was  to  test  the  use  of  cavities  to  decouple  underground  explosions  [36,  42,  43,  44], 

The  subject  of  decoupling  itself  Is  outside  of  the  scope  of  this  report,  but  the  results  of  some 
of  the  Investigations  carried  out  In  connection  with  COWBOY  are  directly  related  to  the  topic. 

Data  obtained  from  some  of  the  shots  In  this  series  apply  directly  to  this  discussion;  the 
shots  were  fired  in  cavities  of  sufficient  size  that  the  response  of  the  salt  was  nearly  elastic, 

If  not  exactly  so.  The  conditions  of  the  theory  outlined  above  were  therefore  satisfied  at  the 
cavity  boundary.  Comparisons  of  theoretical  waveforms  (based  on  clastic  response)  with 
observed  ones  are  presented  by  Parkin  [27]  anc>y  Herbst,  Werth,  and  Springer  [43], 

Parkin  used  the  theoretical  pressures  on  the  cavity  boundaries  calculated  by  Brode  [26]  as 
the  input  to  his  computation  of  displacement  and  other  motion  parameters  as  functions  of  range 
and  time.  The  results  for  the  decoupled  shots  were  excellent;  those  for  the  "overdriven"  cavi¬ 
ties  showed  departures  which  could  be  ascribed  to  plastic  deformation.  Parkin  obtained  his 
results  by  the  direct  numerical  evaluation  of  the  integral  solution  lor  the  displacement  for  an 
arbitrary  pressure  waveform  on  the  cavity  wall. 

Herbst  et  al,  [43]  present  a  theoretical  waveform  for  one  of  the  cavity  shots  In  the  COWBOY 
series  which  agrees  very  closely  with  that  observed.  The  technique  used  here  also  was  numer¬ 
ical  Integration  of  the  equations  for  elastic  response. 

Ah  shown  earlier  In  this  report,  the  smallest  radial  distance  at  which  the  stress  is  low 
enough  for  the  response  to  be  elastic  Is  a  function  of  the  medium,  the  depth  of  the  source,  and 
the  yield.  For  a  given  medium  and  depth,  the  peak  pressuye,  p^^,  at  the  boundary  between  the 
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elastic  region  and  the  nonlinear  zone  is  a  constant.  The  basis  of  decoupling  a  shot  is  that  the 
radius  at  which  the  pressure  is  down  to  Is  much  smaller  for  a  cavity  shot  than  for  one  tightly 
tamped  in  the  medium  [44],  and  therefore  the  radiated  ainplitudo  is  found  to  be  smaller  from 
solutions  like  Equation  40. 

6.2.  THE  SPECTRAL  CONTENT  OF  EXPLOSrON-GENERATi;0  P  WAVES 

Solutions  of  the  form  of  Equations  39  or  40  suggest  that  a  natural  frequency  determined 
by  the  properties  of  the  medium  and  the  size  of  the  charge,  will  characterize  the  compressiv/iia! 
wave.  However,  this  conclusion  ignores  the  fact  that  the  earth  acts  as  a  low-pass  filter  upon 
elastic  waves  [29,  44].  Using  numbers  suggested  by  O'Brien  [7]  for  a  ISO  lb  shot  in  a  medium 
with  a  P-wave  velocity  of  5000  fps  and  Poisson's  ratio  of  0.25,  we  find  that  is  471  rad/sec, 
corresponding  to  a  frequency  of  75  cps.  Fre^ouencies  in  this  range  are  commonly  obser'/ed  in 
the  reflection  method  of  seismic  prospecting,  but  they  are  not  encountered  tn  observations  at 
larger  distances.  Even  if  they  are  transmitted,  most  of  the  instrumentation  used  for  recording 
at  distances  of  more  than  a  few  km  will  not  pass  frequencies  higher  than  about  10  cps. 

Therefore  it  is  necessary  to  examine  the  amplitude  spectrum  corresponding  to  the  solution 
above,  and  attempt  to  determine  the  effect  on  the  observed  amplitude  of  the  low-pass  character 
of  the  transmission -recording  system.  This  has  been  done  by  several  investigators,  all  of  whom 
arrive  at  similar  conclusions  [7,  24,  45,  46,  47]. 

Peet  has  developed  a  theory  In  his  attempt  to  relate  amplitude  to  yield,  in  which  the  shock- 
wave  and  elastic -wave  regions  are  in  Juxtaposition.  The  validity  of  his  results  for  the  elastic- 
wave  region  is  not  in  any  way  affected  by  the  fact  the  terminal  nonlinear  region  has  been  neg¬ 
lected,  «is  long  as  the  actual  radial  distance  to  the  elasvtc  zone  increases  as  the  cube  root  of  the 
yield.  Only  the  numerical  value  radius  of  the  equivalent  cavity  in  terms  of  the  yield  may  be 
Incorrect.  In  turn,  this  affects  only  the  numerical  values  of  amplitude. 

On  the  basis  of  the  output  from  the  hydrodynamic  region,  Peet  finds  the  following  relations 
at  the  input  to  the  elastic  wave  region,  r  =  a: 

p(a,  t)  =  Pj^  exp  ^-Cjt/W^^^  (42) 

where  a  » 

p^  ■  radial  stress  at  which  elastic  behavior  begins  (a  property  of  the  medium) 
t  ■  time  from  start  of  shock  wave 
Cf,  C2  ■  constants 

n  »  empirical  quantity  » 1.16 
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The  pressure  given  by  Equation  42  Is  a  very  much  simplified  verolon  of  all  that  was  covered  in 
the  first  two  sections  of  this  report. 

Because  the  stress  pulse  maintain:  its  siiape  as  it  piopagates  for  pel  '''ct  rla  jtlclty  while 
the  various  kinematic  quan..tie«  particle  velocity)  do  not,  Peet  works  with  the  spectrum 
of  the  mean  stress. 

If  the  displacement  potential,  such  that  u(r,  t)  =  — ,  is  taken  to  be  <>  -  (i  -  the  mean 

stress  ("pressure")  is  p  *  f  ”  t  For  the  pressure  function  of  Equation  42,  it  follows  (for 

1/3 

a  linear  elastic  medium)  that  the  stress  scales  in  both  time  and  amplitude  as  W  '  :  p(r,  t)  = 

.  The  scaling  of  the  pressure  waveform  is  shown  in  Figure  21,  on  which  Q  =  V.h 

The  behavior  of  the  amplitude  spectrum  [24,  Equation  A. 11.13]  le  given  in  Figure  22.  The 

-1/3 

spectrum  has  a  single  maximum  at  a  frequency  f  =  cW  '  ,  at  which  the  amplitude  is  propor¬ 
tional  to  (The  negative  sign  In  the  exponent  is  omitted  from  Peet's  Equations  A.ll.12  and 

A.n.21,  but  should  be  there  from  his  Equation  A.n.20;  see  Figure  23.)  If  the  radius  "a"  increases 

1/3  2 

as  W  '  (as  is  required  by  similitude  considerations  [46]),  the  spectral  maximum  varies  as  a  , 

in  agreement  with  the  pulse  amplitude  in  Equation  40.  For  low  frequencies,  the  spectral  ampli¬ 
tude  is  proportional  to  while  at  very  high  frequencies,  the  amplitude  varies  as  f 

The  effect  of  charge  size  on  the  shape  of  the  spectrum  is  demonstrated  by  Figure  23.  Small 
charges  give  a  flat  apwctrum  with  the  maximum  at  a  high  frequency.  The  difficulty  in  comparing 
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FIGURE  22.  GENERAL  PU-i  Of  1  HE  FREQUENCT  SPECTRUM  OF  pfr,  t) 

/After  Peet) 


FIGURE  23.  THE  FREQUENCY  SPECTRUM  OF  p(r,  t)  FOR 
SEVERAL  CHARGES  (After  Peet) 


amplitude-yield  data  for  amall  and  larf/e  yields,  brought  on  by  the  low-pass  character  of  Ute 

earth -instrument  system,  is  Illustrated  in  Figure  24.  For  small  charges,  the  portion  of  the 

spectrum  which  would  be  recorded  Is  that  for  which  the  dependence  on  yield  Is  near  the  first 
4/3  2/3 

power;  e.g,  from  W  to  W  '  .  For  the  large  charge,  however,  the  spectrum  falls  almost 
entirely  within  the  pass  band,  and  the  exponent  of  the  yield  will  be  smaller. 
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FIGURE  24.  THE  POSITION  OF  THE  FREQUENCY  SPECTRA 
OF  BIG  AND  SMALL  CHARGES  WITH  RESPE'JT  TO  THE 
FASSBANiS  OF  A  LOW-PASS  FILTER  (After  Peet) 

Werth  and  Herbst  [48]  have  ^und  a  first-power  dependence  of  head-wave  amplitude  on 
yield  for  low  yields  and  a  two-thirds  power  dependence  for  high  yields.  The  yield  at  which  the 
dependence  changes  is  higher  for  hard  media  (salt  and  granite)  than  for  soft  ones  (alluvium  and 
tuff). 

Latter,  Martlnelli,  and  Teller  [46]  have  also  analyzed  the  low-frequency  part  of  the  spectrum, 
and  concluded  that  the  signal  at  a  distance  will  vary  as  the  first  or  4/3  power  of  the  yield  de¬ 
pending  on  whether  the  pressure  at  the  elastic-inelastic  boundary  decays  slowly  or  rapidly 
compared  to  the  cutoff  frequency  of  the  earoi-lr.strulnent  system.  Data  precise  enough  to  per¬ 
mit  these  two  scaling  laws  to  be  dlstinguisheu  >i’e  difficult  to  obtain. 

Rather  than  work  from  the  spectrum  for  the  exponentially  decaying  pressure  wave,  O'Brien 
[7]  uses  Blake's  result  [5]  for  radiation  of  a  steady-state  sinusoid  from  the  cavity.  He  finds  that 
the  criterion  for  determining  the  type  of  yield  dependence  Is  whether  wa/c  Is  less  than  or  greater 
than  1.0.  For  low  frequencies,  the  signal  amplitude  (particle  velocity)  should  vary  directly  as 
the  charge  weight.  However,  for  very  large  charges  (a  Is  large),  or  for  very  high  frequencies, 
the  signal  amplitude  will  Increase  as  the  cube  root  of  the  yield.  For  the  first  case,  all  frequen¬ 
cies  will  Increase  in  the  same  manner,  so  that  the  pulse  shape  does  not  change  with  charge  size. 
For  the  second,  the  pulse  will  lengthen  somewhat  with  charge  size. 

Aspolntedout  by  Herbst  etal.[43]and  Carpenter  etal. [45], O'Brien  nssuniesthat  the  radius  of 
the  equt  valent  cavity  increases  as  the  cube  root  of  the  charge  weight,  but  that  this  cavity  is  determined 
by  a  limiting  condition  for  the  mediuniof  Impulse  per  unit  area  equal  to  a  fixed  value.  Since  impulse 
scales  as  W  '  (see  appendix),  these  assumptions  are  inconsistent.  However,  as  shown  by 
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Carpenter  et  al.,  the  step  function  approximation  for  the  stress  at  the  boundary  also  gives  a 
first-power  dependence  on  yield  for  low  frequencies. 

Proceeding  from  very  general  considerations  oi  the  relation  between  pulse  liud  the 
slope  of  the  spectrum,  Weston  [47  J  points  out  that  for  most  disturbances,  the  total  Impulse  at 
distances  more  than  a  few  wavelengths  is  close  to  zero,  and  the  amplitude  spectrum  of  the  pres¬ 
sure  at  low  frequencies  varies  as  f.  Under  conditions  that  length-  and  time-^s^ale  as  he 
then  shows  that  the  low-frequency  amplitude  is  proportional  to  W.  The  variation  of  ttie  low- 
frequency  amplitude  spectrum  with  the  first  power  of  the  frequency  and  yield  results  from  the 
particular  pulse  shape,  ore  for  which  Jpdt  =  0  at  i  iinlte  time,  but  Jfpdtdt  does  not.  In  the 

case  of  a  p<ilse  for  which  tlie  integral  of  the  impulse  is  also  zero  at  intinity,  the  low-frequencv 

2  4/3 

spectral  amplitude  varies  as  f  and  W  '  ,  as  in  Peet’s  result.  Other  exponent.^  for  f  rnd  W  are 
appropriate  for  other  pulse  shapes. 

In  applying  these  general  results  to  underground  explosions,  Weston  looks  for  the  type  of 
source  conditions  that  would  be  consistent  with  these  relations,  for  the  equivalent  cavity  model 
of  an  explosion.  A  step  in  pressure  is  the  simplest  loading  that  will  give  the  prescribed  first- 
power  dependence  of  amplitude  on  frequency  at  a  distance.  The  actual  decay  in  pressure  in  the 
cavity  will  not  affect  the  results  as  long  as  it  is  slow  compared  to  the  longest  waves  of  interest. 
For  a  fast  pressure  decay,  spherical  elastic-wave  theory  may  not  apply  right  up  to  the  radius 
at  which  elastic  behavior  begins,  and  the  cavity  model  may  have  to  be  abandoned.  Weston  states, 
in  conclusion,  that  the  first-power  dependence  Of  amplitude  on  yield  applies  for  frequencies  much 
less  than  1000  cps,  and  for  ranges  greater  thaii  about  700  W^^^  ft  (W  in  lb). 

Carpenter,  Savill,  and  Wright  [45]  have  recently  examined  the  amplitude -yield  relation  in 
1/3 

terms  of  W  '  scaling.  After  considering  the  results  for  a  single  frequency  (which  would 
correspond  to  the  response  of  a  very  sharply  tuned  seismograph),  they  conclude  that  the  band¬ 
width  of  the  instrument  must  be  taken  into  account  when  determining  this  relation  from  field 
data.  As  they  point  out,  and  as  was  seen  above,  different  results  can  be  obtained  for  different 
passbands  of  the  earth-instrument  system.  They  conclude  that  a  simple  W^'^  dependence  is 
likely  to  apply  to  any  particular  frequency  band  up  to  some  limiting  value  of  W,  beyond  which  the 
increase  is  less  rapid. 

These  authors  also  point  out  an  inherent  difficulty  in  applying  any  kind  of  scaling  laws  to 
wave  generation  and  propagation  in  the  real  earth.  Scaling  is  assumed  for  wave  generation;  for 
transmission,  however,  a  constant  geometry  is  assumed,  and  these  are  incompatible  except  for 
an  Infinite  homogeneous  medium.  The  opinion  hold  by  Carpenter  et  at.  Is  that  this  conflict  Is 
not  of  primary  Importance  In  many  practical  situations. 
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6.3.  ENERGY  AVAILABLE  FOR  ELASTIC  WAVES 

The  fraction  of  energy  originally  released  that  is  available  at  the  beg<rning  of  the  elastic 
zone  can  be  calculated  by  the  methods  outlined  for  the  high-pressure  retiions,  onre  adequate 
data  on  the  pertinent  propertl*'^  of  the  medium  have  been  collected.  The  actual  energy  that  Is 
transmitted  to  long  distances  can  only  be  determined  by  Inriud'ng  effects  of  absorption  and 
scattering. 

Parkin  [27]  presents  a  table  of  c.ilculated  values  of  released  and  radiated  e'.<ergy  for  TNT 
shots  of  20  lb  to  2000  lb  fired  in  various  sized  cavities  in  salt.  The  values  are  based  on  Brode's 
results  [26].  The  fraction  of  the  energy  radiated  varies  from  4  x  lo"®  (for  a  completely  de¬ 
coupled  shot)  to  8  X  lO"^  (for  an  overdriven  cavity).  No  estimates  for  tamped  shots  are  given. 

Howell  and  Budenslein  [49]  performed  energy  calculations  for  seismic  waves  generated  by 
explosive  charges  of  about  1.0  lb  at  depths  of  about  10  ft  In  nncnnsohdated  material  tamped  with 
water.  Their  result  for  the  total  energv  in  the  wave  train  (not  the  P  wave  alone)  for  a  surface 
distance  of  10  ft  indicates  that  about  5%  of  the  source  energy  is  present  in  the  frequency  range 
of  six  to  120  cps.  O'Brien  [7]  estimates  that  the  energy  in  the  low  frequencies  Ls  no  more  than 
3%  of  the  source  energy  for  an  underground  explosion,  and  may  be  much  less.  In  view  of  the 
fact  that  the  low-frequency  amplitude  spectrum  Increases  as  the  first  or  second  power  of  the 
frequency  (energy  as  the  second  or  fourth  power),  the  energy  reaching  large  distances  is  likely 
to  be  much  less  than  1%  of  the  yield. ;  , 

Berg  and  Cook  [50]  calculated  the  total  energy  In  the  wave  trains  recorded  at  distances  of 
6.78  to  22  km  from  two  very  large  quarry  blasts;  by  extrapolation  back  to  the  source,  they 
estimated  that  2.7%  of  the  total  energy  content  ^f  the  explosive  charge  appeared  in  the  form  of 
elastic  waves. 

Nlcholls  [28]  used  the  known  properties  of  the  explosive  material  to  determine  yield 
energy,  and  calculated  the  energy  transmitted  into  the  medium  from  strain  measurements. 

His  calculations  indicate  that  2%  to  4%  of  the  chemical  energy  is  coupled  to  the  rock,  but  he 
believes  that  values  of  10%  to  20%  may  be  more  realistic.  The  measurements  were  made  from 
experiments  using  three  different  types  of  dynamite  and  granular  TNT,  tightly  tamped  in  salt. 

He  quotes  the  results  of  earlier  Bureau  of  Mines  investigations  for  granite,  in  which  10%  to  18% 
of  the  explosive  energy  is  coupled  to  the  earth.  These  values  seem  far  higher  than  any  others 
reported. 

Working  with  the  RAINIER  event,  Porzel  [19]  finds  that  only  3%  of  the  original  energy  was 
left  at  the  termination  of  the  crushing  stage  (130  ft),  and  concludes  that  roughly  09%  of  the  energy 
is  deposited  as  heat  in  the  cavity  knd  in  a  thin  shell  of  compressed  material  near  the  source. 
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These  figures  are  In  general  agreement  with  estimates  given  above.  Porzel  concludes  also  that 

the  energy  available  to  produce  seismic  waves  "bears  no  direct  relation  with  the  original  energy 

release,"  but  is  determined  by  the  geometry  of  the  cL.trge  and  hole  in  which  ihe  charge  is  placed, 

and  the  properties  of  the  medium  ( wiporization  pressure  and  dynamic  crushing  jUength).  At 
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long  distances  the  "apparent  energy"  wlil  be  of  the  order  of  10  of  the  original  yield.  This 
value  agrees  with  that  cited  by  Griggs  and  Press  [51]  for  a  20-kt  underground  nuclear  explosion 
(presumably  the  BLANCA  event). 

Although  the  statements  above  are  In  full  accord  with  the  concepts  that  have  been  developed 
In  this  report,  Porzel  arrives  at  some  additional  conclusions  that  are  somewhat  startling.  He 
doubts  that  the  acoustic  approximation  ever  applies  (in  the  sense  of  constant  energy  in  the  wave) 
and  that  it  is  a  misnomer  to  talk  about  partition  of  energy  between  the  source  region  and  the 
radiated  signal.  The  available  energy  decreases  constantly  with  distance  and  rapidly  approaches 
zero,  because  of  known  departures  front  homogeneity  and  linear  elasticity  .  His  conclusion  is 
that  even  the  1%  of  the  explosion  energy  available  at  the  termination  of  the  nonlinear  zone  be¬ 
came  insignificant.  He  suggests  that  the  seismic  signals  "alleged  to  be  observed  as  a  result  of 
the  RAINIER  event"  are  the  result  of  minor  earthquakes  triggered  by  the  passage  of  the  shock 
wave  through  jointed  and  faulted  material.  By  this  reasoning,  the  energy  observed  on  the  seis¬ 
mogram  is  principally  tectonic  strain  energy  accumulated  in  the  test  site  through  geologic  time. 

Press  and  Archambeau  [52]  have  examined  this  hypothesis,  and  present  convincing  argu¬ 
ments  that,  although  an  explosion  will  release  accumulated  strain,  the  magnitude  of  the  energy 
involved  is  too  small  to  affect  the  resulting  seismic  waves. 

Porzel's  argument  that  the  duration  of  the  pi  'ssure  pulse  is  too  short  to  be  resolved  by 
seismic  instruments  overlooks  (1)  mechanisms  in  the  terminal  nonlinear  region  which  act  to 
lengthen  the  pulse  somewhat,  and  (2)  the  fact  that  the  observed  seismogram  will  not  be  the 
pulse  incident  on  the  elastic  wave  region^  but  the  response  of  the  earth -instrument  system  to 
this  pulse.  In  fact,  a  nuclear  explosion,  because  of  its  short  duration,  is  valuable  to  seismolo¬ 
gists  in  that  it  provides  a  direct  observation  of  the  impulse  response  of  the  transmission  path 
and  seismograph,  which  in  turn  can  be  used  to  study  earthquake  focal  mechanisms  [51,  53]. 

The  energy  radiated  in  the  form  of  seismic  waves  depends  upon  many  factors  besides  the 
yield  of  the  explosion;  Parkin's  recent  work  [27]  on  correlation  oi  various  parameters  charac¬ 
terizing  elastic  waves,  therefore,  merits  careful  consideration.  He  correlates  the  theoretical 
peak  values  of  displacement,  particle  velocity,  acceleration,  and  stress  in  salt  with  scaled  range 
by  employing  a  variation  on  established  scaling  procedures.  In  the  following  discussion,  it  is 
assumed  that  the  medium  behaves  elastically. 
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Parkin  first  correlates  the  calculated  values  of  radiated  enerfy  with  the  energy  released  by 
the  explosion.  He  finds  a  good  correlation  of  radiated  energy  with  the  ratio  of  released  energy 
to  cavity  radius.  The  mechanism  of  energy  release  enters  as  a  parameter.  Certainly  a  gen¬ 
eralized  law,  in  which  the  elastic  constants  of  the  medium  are  included,  wc-uld  be  of  great  value. 


Parkin  next  suggests  that  even  though  radiated  and  released  energy  can  be  correlated, 

these  two  quantities  are  of  such  greatly  different  magnitudes  Qat  it  is  worthwhile  to  consider 
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two  different  scale  factors.  For  those  effects  directly  related  to  the  explosion,  W  '  ,  the  cube 

root  of  the  yield  is  to  be  used,  as  In  conventional  procedures;  for  quantities  associated  with 
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propagation  in  the  elastic  medium,  W*  '  ,  the  cube  root  of  the  radiated  energy  is  introduced. 


He  presents  correlations  of  the  various  quantities  mentioned  above.  Only  the  peak  dis¬ 
placement  gave  a  good  correlation  in  dimensionless  form.  By  testing  vai  lous  combinations  of 
dimensionless  velocity,  acceleration,  and  stress  with  ihn  Scale  factor  for  the  elastic  medlUhi, 
he  was  able  to  get  quite  satisfactory  correlations;  the  poorest  was  for  the  acceleration.  He 

finds  that  if  the  dimensionless  quantities  — Jtx,  |r/c,  and  are  multiplied  by 

W  '  c 

n  =  0, 1,  2  for  displacement,  velocity,  or  acceleration,  respectively,  the  resulting  quantities  are 

inversely  proportional  to  the  scaled  range  r/W*^®.  Dimensionless  radial  stress,  ff-VE,  when 

-1/3  1/3 

multiplied  by  W  '  also  is  Inversely  proportional  to  r/W  '  ,  E  is  Young's  modulus. 


This  theory  was  tested  by  using  experimental  values  of  peak  particle  velocity  from  the 
COWBOY  cavity  shots.  The  exponent  of  scaled  range  turned  out  to  be  -1.3,  rather  than  1.0, 
He  suggests  thiit  the  disagreement  can  be  explglnod  by  dissipation  of  energy  due  to  plasticity 
and  other  inelastic  effects  in  the  salt  at  small  ti\rains. 


6.4.  SYMMETRY  OF  THE  SOURCE 

The  entire  development  of  the  theory  of  explosion-generated  waves  has  proceeded  on  the 
assumption  that  the  source  prodiiOfes'a  spherically  symmetrical  stress  wave.  Very  few  obser¬ 
vations  of  the  azimuthal  distribution  of  wave  amplitudes  have  been  published,  and  some  of  these 
are  difficult  to  interpret  because  of  known  variations  in  geological  structure  between  the  source 
and  various  observing  points  and  differences  in  the  geology  at  the  various  recording  sites. 

The  vertical  and  longitudinal  components  of  motion  recorded  at  short  range  from  small 
explosions  in  soil  have  been  found  to  show  negligible  variation  with  azimuth  [10,  31].  On  the 
other  hand,  marked  variations  in  amplitude  are  found  for  all  components  for  similar  experi¬ 
ments  in  rock  [31]. 
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Denoyeretal.  [54]  report  observatlor®  of  asymmetrical  radiation  of  P  and  S  waves  as 

t 

recorded  at  distances  of  190  to  285  km  in  two  dlfferei.t  directions.  They  p«“t*;iate  that  this 
variation  is  caused  at  least  in  part  by  the  fact  that  the  source  w^s  located  Ir.  lov- velon  'ty 
basin  fill  bordered  by  high-veloe'l>  material.  The  two  azimuths  (“150^  apart)  at  which  the 
data  were  available  were  not  oriented  in  any  special  way  with  r'^Siiect  to  the  axes  of  the  basin. 
The  authors  recognized  that  the  differences  in  local  geology  at  ihc  recording  sites  and  the 
presence  of  complicated  structure  between  the  source  and  these  sites  add  to  the  difficulty  of 
evaluating  the  hypothesis. 

Additional  evidence  or.  the  extent  to  which  Uic  source  Is  asymmetrical  is  that  underground 
explosions  always  generate  prominent  SH-  or  Love-type  waves  [9,  37,  38,  52,  55].  If  the  source 
were  completely  symmetric,  horizontally  polarized  shear  waves  could  not  be  generated  in  ?n 
Isotropic,  homogeneous,  or  horizontally  stratified  mediuRi.  Such  motion  can  bo  explained  by 
mode  conversion  at  geologic  or  topographic  irregularities  along  the  transmission  path  and  by 
reflection  of  P  or  SV  waves  from  nonhorizontal  boundaries.  However,  insistence  that  the 
source  is  symmetrical  requires  the  introduction  of  a  variety  of  ^  hoc  assumptions  to  explain 
a  regularly  observed  phenomenon. 

Processes  in  the  nonlinear  region  around  the  explosion  (discussed  in  earlier  parts  of  this 
report)  seem  to  be  an  adequate  source  of  SH-wave  generation,  and  even  nonsymmetrlc  P-  and 
Rayleigh-wave  generation.  Wright  and  Carpenter  [38]  suggest  cracking  in  hard  rock  and  in¬ 
stability  of  the  gas-soil  boundary  in  unconsolidated  material  as  possible  sources.  Kissllnger  et 
al.[37]have  suggested  radial  crack  production  as  the  source  of  SH  motion  that  they  have  observed. 

Apparently  no  one  but  AokI  has  attempted  to  Cv-Vmlate  the  elastic  waves  proceeding  from 
the  nonlinear  region,  and  he  starts  with  the  assumption  of  spherical  symmetry.  Short  has 
observed  that  cracking  in  salt  does  not  occur  equally  in  all  directions,  but  that  large  cracks 
extend  In  a  few  directions  [56].  If  this  is  generally  true,  it  would  seem  difficult  but  not  impos¬ 
sible  to  attempt  an  analysis  of  the  wave  motion  generated  by  cracking.  The  distribution  of 
first-motion  amplitude  of  P  and  S  waves  around  a  growing  crack  is  given  by  Knopoff  and  Gilbert 
[57],  It  Is  necessary  to  estimate  the  number  and  direction  of  cracks,  their  rate  of  growth  in 
terms  of  the  properties  of  the  medium;  and  the  amplitude  of  displacement  normal  to  the  crack. 

* 

IXPIRIMINTAl  STUDIiS  OF  WAVE  OENERATION  BY  EXPLOSIONS 

.  ;  fir.  ■  • 

Knowledge  concerning  the  characteristics  of  explosion-generated  seismic  waves  is  based 
primarily  on  observations,  because  of  the  many  uncertainties  Involved  in  the  behavior  of  the 
explosion,  the  coupling  to  the  earth,  arid  the  properties  of  the  medium.  In  preceding  sections 
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some  empirical  relations  have  been  Included  along  with  the  theoretical  development  presented 
In  order  to  Illustrate  particular  points.  This  section  presents  a  brief,  ov^iematlc  review  of 
published  observations. 

It  is  not  practical  to  discuss  ec.'h  published  paper  on  this  subject  In  detail,  owing  to  the 
large  number  In  existence.  Fortunately,  a  number  of  re ’dev  papers  consolidate  information 
from  several  Investigations;  such  reviews  'vill  be  the  primary  source  of  the  material  presented 
herein.  Each  review  paper  contains  an  appropriate  bibliography,  and  the  reader  interested  in 
the  details  of  experimental  methods  and  conditions  will  find  it  advisable  to  go  back  to  the  orig  ¬ 
inal  papers. 

Experimental  studies  are  valuable  for  gaining  insight  Into  the  phenomena  Involved  and  for 
predicting  the  results  of  further  evpffr'.ments.  In  the  absence  of  a  complete  theory,  it  is  im¬ 
portant  to  be  able  to  treat  a  given  experiment  as  a  mode!  of  all  experiments.  This  requires  an 
understanding  of  how  the  effects  of  the  explosion  scale  with  the  various  parameters  that  char¬ 
acterize  the  experiment.  Some  scaling  relations  have  been  presented  In  the  previous  section, 

A  brief  discussion  of  the  principles  of  scaling  is  presented  in  the  appendix. 

The  following  discussion  covers  effects  of  depth  of  burial,  medium,  and  the  combination  of 
yield  and  distance. 

7,1.  THE  EFFECT  OF  SOURCE  bEPTH 

When  the  effect  of  the  depth  of  burial  on  the  resulting  seismic  signal  is  considered,  it  is 
necessary  to  distinguish  between  .effects  on  lx  S’  waves  (the  P  wave  In  particular)  and  effects 
on  surface  waves.  Dy  tlielr  very'nature,  surface  waves  depend  for  their  existence  on  the  pres¬ 
ence  of  a  free  surface,  and  are  sb'ongly  affected  by  the  distance  of  the  source  from  this  sur¬ 
face  [58,  59].  Several  studies  of  ijurface-wave  generation  and  propagation  have  been  published, 
but  they  will  not  be  reviewed  here  [59,  60,  61,  62], 


There  is  no  such  direct  effect  for  body  waves.  It  Is  apparent  from  the  theory  presented 
above,  however,  that  the  radial  distance  at  which  elastic  behavior  begins  Is  a  function  of  the 
llthostatlc  (or  ambient)  pressure^  In  turn,  this  radius  has  been  shown  to  directly  affect  the 
signal  amplitude.  Bigger  "equlv^ent  cavities"  give  bigger  signals.  In  a  homogeneous  medium, 
the  P-wavo  amplitude  can  be  expected  to  Increase  rather  sharply  as  the  depth  of  the  charge  Is 

■t'i':  .  -V 

increased  from  zero  to  that  required  for  complete  containment,  and  then  decrease  gradually 
for  greater  shot  depths  (lncreas<l^  atnblcnt  stress).  At  most  test  sites.  It  Is  very  likely  that 
beds  of  rock  of  differing  propert^s  arc  encountered  as  the  shot  depth  Is  Increased,  so  that 
depth  alone  Is  not  the  controlUnB^nqtor. 
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Experiments  on  the  effect  of  source  depth  on  the  P  wave  are  subject  to  soit.a  Inherent  dif¬ 
ficulties.  If  the  observations  are  made  ai  short  ranges  with  fixed  detector  positions,  the  effects 
of  source-detector  distance  may  overwhelm  the  effect  of  depth.  Since  absorption  losses  are 
Imperfectly  known,  corrections  for  v-oMatlons  In  distance  are  uncertain.  If  a  se^srauineter  Is 
located  on  the  surface,  another  problem  arises.  As  Is  well  known,  the  surface  displacement 
Is  not  the  same  as  the  Incident  displacement,  but  is  a  resubaal  of  the  incident  and  reflected 
waves.  The  observed  displacement  Is  a  complicated  function  of  the  angle  of  incidence  and  the 
elastic  properties  of  the  medium  [63J.  At  short  ranges  for  which  the  angle  of  Incidence  d**  - 
creases  rapidly  with  Increasing  depth,  therefore,  the  data  must  be  Interpreted  with  care.  The 
use  of  deeply  buried  detectors  eliminates  this  pr.'blem. 

On  the  other  band,  when  observations  are  made  at  distances  great  enough  so  that  the  ancle 
of  Incidence  Is  not  greatly  affected  by  bhe  rU>pth,  normal  geologic  conditions  will  usually  result 
In  the  first  arrival  being  a  head  wave  from  some  subsurface  boundary.  O'  course,  the  effect  of 
depth  on  this  head  wave  may  be  of  Interest,  but  examination  of  the  direct  P  wave  becomes  dif¬ 
ficult. 

Stengel  has  studied  shot-depth  effect  for  1. 5-lb  charges  of  dynamite  [64];  the  principal  results 
are  reported  In  a  more  easily  accessible  paper  by  Howell  [65].  The  depth  was  varied  from  7 
to  102  ft,  and  the  media  were  dolomite,  weathered  dolomite,  and  gravel  and  sand  near  the  sur¬ 
face.  A  close-in  Instrument  was  moved  as  the  depth  changed  to  keep  the  angle  of  Incidence 
roughly  the  same,  and  an  Instrument  at  200  ft  from  the  shot  hole  was  kept  fixed.  Relative 
energy  In  the  wave  train  was  calculated  from  particle  velocity  records. 

The  energy-depth  data  for  the  close  seismon  ater  oscillate  through  an  order  of  magnitude 
for  changes  In  depth  of  about  10  ft.  The  extreme  excursions  are  bounded  by  parallel  trends 
corresponding  to  an  exponential  increase  of  energy  with  depth,  given  by  E^  »  E^  exp  (0-047  d). 
The  total  relative  energy  at  200  ft  Is  described  quite  well  by  the  same  relation.  For  the  deeper 
shots  a  strong  refraction  arrival  Increased  In  energy  at  a  slightly  higher  rate.  Obviously,  an 
exponential  Increase  of  relative  energy  with  depth  can  be  valid  only  for  relatively  shallow 
depths. 

In  his  pioneering  study  of  blast  effects  In  soli,  Lampson  [30]  found  that  at  scaled  ranges  of 
r/W^/3  between  2  and  16,  both  pressure  and  positive  Impulse  per  unit  area  were  linearly  de¬ 
pendent  on  a  coupling  factor  wtiich  depended  on  depth.  This  factor  for  one  of  the  soils  In  which 
he  worked  Is  presented  In  Figure  25.  The  factor  has  a  maximum  at  a  depth  of  about  2W^'^^ 

(W  In  pounds):  It  falls  off  rapidly  for  shallower  depths,  and  more  slowly  for  great  depths. 
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FIGURE  25.  EXPLOSIVE  COUPLING  FACTOR  AS  A  FUNCTION  OF  DEPTH  OF  CHARGE: 
CLAYEY  SILT  SOIL  (After  Lampson) 


Rocard  [66]  has  recently  published  results  of  a  preliminary  attempt  to  combine  data  from 
a  v-dde  variety  of  sources  Into  a  single  curve  pf  relative  seismic  amplitude  at.  long  ranges  as 
a  function  of  depth  (see  Figure  26).  The  data  points  b  are  from  numerous  small  charges  re¬ 
corded  at  5  km;  n  indicates  nuclear  explosions  of  the  HARDTACK  II  series;  m  represents 
underwater  explosions  recorded  at  110  and  300  km;  and  g  means  submarine  grenades  recorded 
at  160  to  190  km.  All  data  were  scaled  to  a  charge  wplght  of  1.0  kg  by  dividing  the  depth  by 
the  cube  root  of  the  charge  we,lght  In  kilograms.  The  curve  was  developed  from  pairs  of  data 
points  corresponding  to  similar  recording  conditions  and  different  scaled  depths. 

‘  h* 

The  unit  of  "seismic  effect"  bn  this  curve  was  taken  as  the  signal  from  a  half -burled 
charge.  Rocard  suggests  (private  communication)  th4t  an  e.-^ploslon  that  Is  just  deep  enough 
to  be  contained  Is  a  better  experimental  standard.  Such  an  explosion  has  a  value  of  about  5  on 
the  curve,  so  the  abscissas  should  be  divided  by  5  to  give  the  relative  seismic  effect  In  terms 
of  this  revised  unit.  The  curve  shows  an  Increase  by  a  factor  of  ~6  at  great  depth  relative  to 
the  barely  contained  charge. 
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Altitude 


FIGURE  26.  RELATIVE  AMPLITUDE  OF  SEISMIC  SIGNAL 
AS  A  FUNCTION  OF  DEPTH  OR  ALTITUDE  (After  Rocard) 


The  vertical  curve  for  surface  and  above-giuund  detonations  of  small  charges  and  the 
dtscoTsttnuity  at  depths  Just  below  the  surface  agr  ^  roughly,  at  least,  with  recent  unpublished 
observations  of  the  present  author. 

Ricker  [67]  found  an  Increase  in  the  seismic  efficiency  of  buried  charges  to  a  depth  of 
about  100  ft  for  1.0-lb  charges  of  dynamite  from  his  experiments  in  the  Pierre  shale.  Scaled 
to  1.0  kg,  this  maximum  is  at  a  depth  of  40  meters,  a  result  not  discordant  with  Rocard's 
curve. 

Rocard  has  pointed  out  (private  communication)  that  the  data  on  the  curve  (Figure  26)  for 
the  deeper  shots  were  obtained  from  tests  in  the  Mediterranean  Sea  that  were  deep  compared 
to  the  barely  contained  expolosion,  but  well  above  the  sea  bottom.  No  data  are  Included  from 
shots  so  near  the  bottom  that  interference  from  that  boundary  might  alter  the  results. 

Adams  and  Swift  [35]  found,  for  two  TNT  shots  in  tuff,  that  the  radiated  elastic  energy 
tor  a  shot  nt  a  depth  of  1010  ft  was  roughly  four  times  that  for  a  shot  at  240  ft.  The  measured 
ambient  stress  was  five  times  as  great  for  the  deeper  explosions.  When  scaled  to  1.0  kg, 
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these  depths  are  44  and  9.4  meters,  and  an  amplitude  Increase  by  a  factor  ot  2  is  concordant 
with  an  extension  of  Rocard's  curve. 

7.2.  THE  EFFECT  OF  THE  MEDH’M 

The  manner  In  which  the  properties  of  the  medium  enter  Int..  'he  wave-generation  process 
was  indicated  in  the  sections  concerning  theoretical  development.  The  Hugonlot  equation  of 
state,  the  crushing  strength,  elastic  limit,  porosity,  fluid  content,  and  brltt’rness  all  play  roles 
Implied  by  the  equations  describing  the  progress  of  the  pressure  pulse. 

The  value  of  Poisson's  ratio  plays  a  key  role  In  the  elastic-wave  region  at  short  range, 
determining  whether  the  pulse  is  highly  oscillatory  cr  pulse-like.  The  slxe  of  the  inelastic 
region  and  tlic  peak  pressure  at  Its  boundary  (both  dependent  on  the  medium)  deterriine  tl.t 
amplitude  for  a  given  yield.  The  predointnsnt  period  in  the  P  wave  dopends  on  tht  size  of  the 
equivalent  cavity,  the  compreRstQr.al  «ave  velocity,  and  Polssun’s  ratio. 

As  suggested  by  the  quotation  from  Sharpe  at  the  beginning  of  this  section,  a  shot  in  weath¬ 
ered  material  produces  a  broader  pulse  than  one  in  unweathered  material,  and  ?  shot  in  hard 
rock  produces  a  sharper  pulse  than  one  in  softer  material.  This  has  been  confirmed  experi¬ 
mentally  [68].  Field  records  illustrating  the  effect  of  the  medium  at  short  ranges  are  shown 
In  Figure  27  (69).  The  limestone  records  are  characterized  by  much  higher  frequencies  than 
the  other  media.  At  distances  over  100  meters,  the  P  wave  is  attenuated  below  the  level  of 
sensitivity  of  the  instrumentation  and  the  records  are  dominated  by  surface  waves,  which  are 
of  much  greater  amplitude  in  the  soil  than  in  rock. 

Although  every  experimental  study  of  expl  «!cn-generated  seismic  waves  includes  the 
effect  of  the  medium,  only  a  few  have  been  devoted  specifically  to  this  problem.  Lampson 
[30,  p.  26]  concludes  that  the  type  of  soil  is  the  most  Important  single  variable  governing  the 
transmission  of  the  stress  wave  from  an  underground  explosion.  His  work  pertains  only  to 
■oil,  of  which  five  different  types  were  Included.  He  defines  a  soil  constant  (with  dimensions 
of  an  elastic  modulus)  which  was  found  to  vary  from  an  average  value  of  800  for  loess  to 
100,000  for  saturated  clay.  The  pressure  at  a  given  scaled  range  was  found  to  depend  linearly 
on  this  constant.  This  "constant"  is  variable  over  a  rather  wide  range  for  a  given  soil,  depend¬ 
ing  on  the  moisture  content  and  the  compaction.  A  reasonably  good  correlation  of  the  soil  con¬ 
stant  and  the  eompressional  seismic-wave  velocity  was  found.  This  is  a  bit  unexpected,  since 
the  seismic  wave  velocity  pertains  to  small  strains  and  the  soil  constant  the  entire  (nonlinear) 
■tress -strain  curve  up  to  the  peak  stress, 

A  recent  set  of  experiments,  specifically  designed  to  test  the  effect  of  the  medium  on 
coupling,  was  carried  out  for  the  two  media  which  have  received  the  greatest  attention  in  this 
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report:  tuff  and  salt  [35).  Charges  o'  grar.ular  TNT  of  100  to  1000  lb  were  tamped  In  holes  at 
depths  such  that  the  ambient  stress  levels  we**e  similar.  Only  the  rel"“.Ta  strengths  of  the 
signals  In  the  two  media  are  Investigated,  with  no  attempt  to  get  an  absolute  fotinf  ^g  coefficient 
for  either.  The  result,  sh'V.r'ng  the  frequency  dependence  of  the  medium  effect,  Is  presented  In 
Figure  28.  The  authors  (Adams  and  Swift)  conclude  that  tft^  ratio  of  signal  strength  In  tuff  to 
that  in  salt  is  1.6  ±  0.4  for  the  frequency  range  of  30  to  60  cpa.  This  corresponds  to  a  ratio 
of  energies  of  roughly  2.6  to  1. 
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FIGURE  28.  MEDIUM  COUPLING  FOR  TUFF  TO  SALT  AS 
A  FUNCTION  OF  FREQUENCY  (After  Adams  and  Swift) 


It  should  be  noted  that  this  result  Is  at  least  partly  explained  by  Nicholl’s  findings  for  the 
effect  of  the  ratio  of  Impedance  of  the  explosive  to  that  of  the  rock  [28].  Nlcholl's  gives  a 
value  for  the  characteristic  Impedance  of  granular  TNT  of  16.7  Ib-sec/ln^  and  for  salt  of  35 
Ib-sec/ln^.  From  the  density  and  compresslonal  velocity  lor  tuff  at  a  depth  of  1000  ft  from 
Adams  and  Swift,  the  Impedance  of  the  tuff  Is  16.6  ib-sec/ln'"'.  Thuu,  the  Impedance  ratios  are 
0.48  for  salt  and  0,99  for  tuff.  Applying  these  values  to  the  average  curve  In  Figure  11  of 
Nlcholl's  paper,  a  ratio  of  energy  coupled  Into  the  tuff  to  that  In  the  salt  Is  about  1.4,  or  a  tuff- 
to-salt  signal  amplitude  ratio  of  1.2. 

Though  the  result  given  by  Adams  and  Swift  Is  a  time  medium  effect,  it  Is  valid  only  for 
the  particular  explosive.  It  even  seems  possible  to  reverse  the  result  by  using  hlph-veloclty 
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dynamite  (Impedance  of  32.3  Ib-sec/tn^),  although  Nicholls  gives  no  discussion  of  behavior  for 
Impedance  ratios  greater  than  1.0.  It  would  seem  ^.'t  future  experiments  on  the  effect  of  the 
medium,  using  chemical  charges,  should  employ  several  types  of  explosivps. 

Werth  and  Herbst  [48]  have  ;.-tudled  amplitudes  of  seismic  waves  generated  by  nuclear  ex¬ 
plosions  In  tuff,  salt,  granite,  and  alluvium.  Ttiey  Isolate  tiie  effect  of  the  medium  at  the  source 
by  including  in  their  analysis  the  effects  of  couisllng  out  of  the  sour4.c  region  and  of  the  crustal 

i 

structure  between  the  source  and  the  recording  site.  The  theory  developed  In  Reference  1  is 
applied.  All  results  are  scaled  to  a  yield  of  5  }t$.  by  the  wl/3  relationship.  The  amplitudes, 
at  distances  of  the  order  of  500  km  relative  to  b,»rf,  are:  biff,  1;  alluvium,  0.25;  granite,  1.11; 
and  salt,  1.61.  If  realistic  models  of  the  crustal  structure  are  inctuded,  the  amplitude  ratios 
for  the  first  half-cycle  of  the  head  wave  from  t^e  Mohorovtclc  discontinuity  are:  tuff,  1; 
alluvium,  0.18;  granite,  2.26;  and  salt.  ?.tl.  The  result  oi[  Adame  and  S  dft  for  the  tuff-to-salt 
ratio  is  seen  to  be  reversed.  '1'  v 

Although  the  experimental  work  under  discussion  Is  ||irlmarlly  field  work,  worthwhile 
Information  on  the  wave-generating  process  can  be  gatnet|from  laboratory  experiments.  Such 
tests  lend  themselves  especially  well  to  work  on  the  properties  of  the  medium,  since  there  is 
an  obvious  limit  to  the  range  of  yields  and  ejtfectlve  source  depths  that  is  practical  to  Investl- 
gate.  ■ 

-I' 

The  basic  principles  of  shock-wave  propagation,  as  outlined  in  the  discussion  of  the 
Rankine-Hugontot  equations,  have  been  used* directly  by  Hughes  and  McQueen  [70]  to  measure 
rock  densities  at  pressures  from  150  to  750  kilubiti.s.  Although  their  objective  was  to  obtain 
information  of  direct  importance  to  geophysical  scales  of  the  earth's  Interior,  It  Is  obvious 
from  the  many  comments  on  the  lack  of  data  pertaining  to  rock  properties  at  high  pressures 
that  experiments  of  this  type  are  highly  important  to  the  solution  of  the  present  problem. 

Ito,  Terada,  and  Sukurat  [71]  used  similar  techniques  In  their  Investigation  of  wave  gen¬ 
eration  in  sandstone  and  marble.  Shock  waves  produced  by  blasting  caps  and  300-gram  high- 
explosive  charges  were  used  to  load  the  specimens.  Stress  was  measured  by  shock-wave 
techniques,  and  strain  was  directly  recorded  by  suitable  gauges. 

The  stress  wave  showed  a  decrease  in  rise  time  with  increasing  peak  pressure  for  marble, 
but  not  for  sandstone.  These  conclusions  are  based  on  peak  pressures  from  blasting  caps  of 
about  0.3  to  1.6  ktlobara.  Dynamic  values  of  Young’s  modulus  were  determined  to  be  two  or 
three  times  greater  than  static  values;' dynamic  st;ength  is  also  greater  than  static  strength. 
Evidence  of  plastic  flow  was  found  for  the  explosive  charges,  but  not  for  the  caps. 
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The  internal  structure  of  the  marMe  specimen  was  altered  much  more  than  that  of  the 
sandstone  for  the  same  loadint;.  This  conclusion  was  reached  by  running  <iiutic.  compression 
tests  on  the  samples  after  the  shock  wave  tests  Direct  evidence  of  the  1:;«’  veiocil,'  of  shock 
propagation  In  certain  parts  o1  Uie  nonlinear  region  was  found. 

The  peak  pressure  was  found  to  decrease  very  rapld'y  *ith  distance  from  the  source,  and 
differences  in  peak  pressures  for  different  kinds  of  explosive  materials  ■sere  evident  only  in 
a  narrow  range  close  to  the  source. 

7.3.  RELATION  OF  SIGNAL  STRENGTH  TO  RANGE  AND  YIELD 

Whether  an  Investigator  is  interested  in  explosion-generated  waves  from  the  viev'polnt  of 
quarry  blasting,  seismic  exploration,  or  detection,  he  obviously  is  concerned  with  the.nanner  in 
which  the  various  parameters  associated  with  the  stress  wave  vary  w  th  the  yield  and  the  dis¬ 
tance  from  the  source.  Consequently,  the  published  studies  on  this  aspect  of  the  problem  are 
numerous,  and  a  complete  review  of  all  the  literature  on  this  subject  will  not  be  undertaken. 

The  effects  of  the  medium  around  the  shot  and  the  depth  of  the  source  enter  into  all  of 
these  studies.  The  results  presented  herein  are  empirical  relations  pertaining  to  the  particu¬ 
lar  conditions  of  the  test  on  which  they  are  based,  and  extrapolation  to  other  situations  must 
be  approached  cautiously. 

A  large  number  of  reports  have  been  published  covering  Individual  projects  connected  with 
the  explosion  program  of  the  U.  S.  Atomic  Energy  Commission.  Most  of  these  are  in  the  form 
of  project  reports  that  are  not  widely  accessible.  However,  papers  reviewing  large  parts  of 
this  work  have  recently  appeared  in  the  regular  rctentific  periodicals,  and  these  will  be  used 
as  the  primary  references.  ■ 

The  work  of  the  U.  S.  Coast  and  Geodetic  Survey  and  several  other  individuals  and  organi¬ 
zations  on  the  measurements  of  surface  motion  is  summarized  by  Carder  and  Mickey  [72]. 

■  i  .  T: 

Surface  displacements  and  accelerations  produced  by  yields  of  180  lb  to  19.2  kllotons  (equiva¬ 
lent)  of  TNT  have  been  measured  at  distances  from  0.1  to  900  km.  These  shots  have  been 
fired  in  salt,  tuff,  alluvium,  limestone,  quartzite,  and  basalt.  It  will  be  noted  that  these  meas¬ 
urements  are  outside  of  the  f^ee-fleld  region,  and  are  Influenced  by  geologic  structure  and  the 
resulting  reflections  and  refractions.  The  results  are  entirely  concerned  with  maximum  values 
of  dtspalcement  and  acceleration,  without  regard  to  the  particular  wave  type  or  path  with  which 
the  values  are  associated.  In  a  sense,  the  results  do  not  directly  relate  to  the  principal  sub¬ 
ject  of  this  report,  the  initial  P  wave.  These  results  are  included,  however,  because  they  are 

t-.  ■. 

certainly  a  measure  of  the  efiocttvqness  of  the  explosion  an  a  seismic  source. 


Peak  acceleration  data  for  all  shots  fired  <n  the  Rainier  Mesa  at  uie  Nevada  Test  Site 
are  summarized  in  Figure  29.  The  ground  acceleration  prediction  formula  is 

a- (2.4  ±1.0)  (43)  ^ 

where  a  =  peak  acceleration  in  units  of  gravity 

W  =  yield  in  kilotons  ' 

r  =  range  In  thousands  of  feet  '''  i*; 

The  data  on  Figure  29  have  been  normalized  to  1.0  kt  by  the  law.  The  formula  is  intended  '  > 

to  be  used  from  outside  the  free  field  zone  to  5  or  6  km.  "/j 


■■ 
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A  general  scaling  law  for  maximum  earth  displacement  has  go;.e  through  a  number  of  re¬ 
visions.  The  form  offered  by  Carder  and  Mickey  is 

A^  CW®  *V‘'x  10‘“  ‘ 

where  A  =  maximum  surface  displacement  (cm) 

C  =  a  constant 

r  =  source-to-detector  distance  (ft) 
n,  k  are  constants 

f  =  dominant  frequency  by  visual  examination  (cps) 

The  values  of  the  constants  have  been  found  to  depend  upon  the  distance  range.  With  reference 
to  the  300  to  flIIOO-ft  range  tn  Figure  30,  C  is  10^/5"^,  n  Is  2,  and  k  is  1,75  x  10"®;  in  the 
9800-ft  to  100-ml  range,  C  is  10®-®,  «  is  1,  and  k  Is  1.83  x  10"®.  From  ll'-O  to  60 J  ml,  C  Is 
Iq-2.82^  n  is  0.5,  and  k  is  7.6  x  10“^.  These  amplitudes  are  multiplied  by  about  3  for  instru¬ 
ments  mounted  on  alluvium. 

The  work  conducted  by  several  agencies  on  strong  motion  from  underground  niirleav  shots 
Is  summarized  in  the  report  by  Adams  et  al.  (2].  Since  these  data  were  taken  In  the  free-fteld 
zone,  they  pertain  directly  to  the  present  report.  The  largest  amplitude  of  acceleration  or 
displacement  tn  the  direct  wave  is  used  tn  the  studies  of  surface  motion. 

Subsurface  peak  radial  acceleration  was  found  to  decrease  as  the  cube  or  fourth  power  of 
the  slant  range.  A  tendency  toward  a  second  or  first  power  dependence  is  Indicated  at  larger 
ranges,  but  Is  masked  by  effects  produced  bv  local  geology.  Adams  et  al.  point  out  that  accel¬ 
eration  data,  especially  those  collected  at  the  v<riace,  are  greatly  affected  by  the  presence  of 
bedding  planes  or  other  surfaces  of  discontinuity  In  the  medium.  They  also  point  out  that  par¬ 
ticle  velocity  may  prove  to  be  a  better  parameter  than  acceleration  to  use  as  a  reference,  be¬ 
cause  It  is  closer  In  waveform  to  the  stress  pulse , 

Vertical,  radial,  and  tangential  components  of  surface  acceleration  were  all  found  to  follow 
the  same  formula  with  satisfactory  precision; 


where  a  ■  units  of  gravity 
W  ■  kllotons 
r  •  slant  range  (ft) 


W®-Vr-2-^ 


X  10 


,8.5 


This  formula  may  be  compared  with  that  of  the  Const  and  Geodetic  Survey  by  expressing  r  In 
thousands  of  ft.  :^he  coefficient  In  Equation  43  becomes  1.6,  near  the  lower  limit  of  the  Coast 
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GROUND  EFFKCtS  FROM  UNDERGROUND  EXPLOSIONS 


riOURE  30,  PEAK  PARTICLE  DISPLACEMENT  VS.  DISTANCE 
NORMA UZED  TO  A  lo”-LB  YIELD  BY  THE  SCALING  FUNC¬ 
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and  Geodetic  Survey  result.  Adams  et  al.  show  that  this  result  Is  iiot  consistent  with  free-fleld 
1/3 

scaling  based  on  W  '  scaling  of  length. 

On  the  other  hand,  the  displacement  data  are  found  to  fit  polnt-sourct  i-  "allrj-'  laws  much 

more  closely.  Each  of  the  l^reo  components  was  found  to  fit  a  relationship  A  =  mw"c"P,  with 

slightly  different  values  for  the  constants.  For  the  vertlca*  component,  n  was  found  to  be  0.69 

and  p  to  be  1.14.  The  corresponding  values  lor  the  radial  coniponent  were  0.94  and  1.68,  If  the 

{n  ^  1^/  3  “T 

exponent  of  r  Is  inserted  Into  the  scaling  relationship  for  displacement,  A  =  kW  r  ,  the 
resulting  exponent  for  W  is  quite  close  to  the  observed  value.  The  authors  feels  that  this  follows, 
tn  spite  of  the  poor  agreement  for  acceleration,  from  the  fact  that  the  displacment  Is  associated 
with  longer  periods  and  is  not  so  greatly  affected  by  the  presence  of  interfaces. 

The  exponent  of  the  yield  is  ch*ea  io  that  of  the  Coast  and  Geodetic  Survey  for  the  vertical 
and  transverse  component,  and  somewhat  closer  to  iitilly  for  the  radial  component.  It  will  be 
recalled  that  an  exponent  of  4/3  or  1.0  was  called  for  at  the  low-frequency  part  of  the  spectrum 
in  the  theoretical  development,  depending  on  the  pulse  shape.  Peet  predicts  a  value  of  0.67  near 
the  maximum  of  the  spectrum.  It  may  be  seen  from  these  results  that  the  observed  relations 
fall  In  the  expected  range,  but  that  the  experimental  data  are  not  precise  enough  to  verify  the 
theory.  This  Is  not  surprising.  In  view  of  the  nature  of  explosions  and  earth  materials. 

Willis  and  Wilson  [73]  have  found  that  the  maximum  vertical  displacement  Increases  as 
the  first  power  of  the  yield.  This  conclusion  is  based  on  observations  of  13  quarry  blasts,  two 
large  excavation  blasts,  and  two  underground  nuclear  explosions.  It  is  well  known  that  correlation 
of  signal  strengths  from  quarry  blasts  are  iilgiiiy  variable  [74],  and  depend  on  such  factors  as 
th9  spacing  and  burden  of  the  shot  and  the  loc;  *von  of  the  shot  in  a  given  quarry. 

On  scaling  all  their  data  to  1.7  kllotons  (RAINIER  yield)  by  means  of  this  W^‘®  relation, 
these  same  authors  have  found  empirical  amplitude-yield  relations  (for  r  between  1  and  200  km) 
as  follows;  A  *  (0,65  ±  0.15)r"^^^,  where  A  maximum  vertical  amplitude  (normalized  to 
1.7  kllotons)  In  cm,  and  r  >  distance  In  km. 

Between  100  and  1000  km,  ' 

A  «  (0.013  ±  0.003)  exp  [-(00.72  ±  0.0003)r]  r*^^^ 

64%  of  the  data  fall  within  the  envelopes  of  these  equations;  the  data  and  empirical  curve  are 
shown  In  Figure  31, 

Additional  attempts  at  evolving  amplltude-yleld-dlstance  equations  are  also  to  l)e  found  in 
Berg  and  Cook  [50],  Oaskell  [75],  Hablmrjam  and  Whotton  [76],  Ito  [77],  Lampson  [30],  Murphey 
[36],  OI)ert  and  Duvall  [78],  O'Brien  [7],  and  Thoenen  and  Wlndes  [79].  Some  of  these  are  ap- 
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FIGimE  31.  PEAK  VERTICAL  GROUND  DISPLACEMENTS  FOR  EXPLOSIVE  BLASTS 

NORMAUZED  TO  THE  RAINIER  EVENT  BY  THE  W*’®  SCALING  FUNCTION  (After 

Willie  end  'vilson) 


plicable  to  specific  media,  e.g.,  Lampson's  work  on  soils  and  Murphey's  in  salt.  Others  are 
primarily  concerned  with  quarry  blasts,  e.g.,  the  pioneering  work  of  Thoenen  and  Windes,  and 
the  studies  of  Habberjam  and  Whetton. 

Fundamental  studies  on  stress-wave  generation  are  described  in  a  long  series  of  papers 
by  investigators  at  the  U.  S.  Bureau  of  Mines.  In  addition  to  the  report  of  Obert  and  Duvall  cited 
above,  work  by  Nlcholls  f28]  and  by  Fogelson  et  al.  [80]  contain  data  on  strain  propagation  in 
rocks.  References  to  earlier  Bureau  of  Mines  studies  tnay  be  found  in  these  papers. 

There  seems  little  value  in  recounting  here  in  detail  the  specific  formulas  at  which  these 
many  investigators  arrived,  since  an  evaluation  of  their  pertinence  requires  a  full  description 
of  the  experimental  procedures  and  conditions  as  well  as  the  methods  of  analysis.  Both  O'Brien 
and  tto  present  summaries  of  some  of  these  studies  in  addition  to  their  original  data.  Yield  dc- 
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pendencies  ranging  from  the  0.4  powe-  to  the  1.1  power  have  been  found.  O'Brien  finds  that, 
for  those  studies  in  which  u>a/c  (see  Section  6.1)  is  known  to  be  less  thor  l.O,  tiie  first  power 
dependency  seems  to  hold. 

Although  empirical  relations  of  signal  strength  and  range  ca'  be  found,  any  study  which 
Ignores  attenuation  resulting  from  Inelastic  behavior  of  medium  is  an  oversimplification. 

8 

SUMMARY  AND  CONCLUSIONS:  THE  STATE  OF  THE  ART 

The  study  of  explosion-generated  seismic  waves  is  stiil  primal  ily  an  empirical  science. 

In  the  absence  of  a  complete  theory  of  the  processes  that  result  in  the  radiated  wave,  progress 
in  further  understanding  of  the  sub)';'}  ‘s  dependent  on  well  designed  experiments.  Uniorlu.uitc- 
ly,  the  complex  properties  of  earth  materials  and  geologic',  structure  in  v.chich  full-scale  tests 
must  be  carried  out  make  the  design  of  experiments  difficult. 

It  has  been  seen  that  the  analytic  tools  needed  to  develop  the  theory  are  known,  and  seem 
to  be  adequate.  Numerical  techniques  exist  for  evaluating  solutions  that  cannot  readily  be  put 
into  closed  form.  The  greatest  need  is  for  additional  experimental  evidence  concerning  the 
nature  of  the  processes  Involved,  so  that  the  theoretical  development  can  be  guided  along  realistic 
lines.  This  Is  particularly  true  for  the  terminal  nonlinear  region.  Although  both  experimental 
and  theoretical  work  in  this  zone  are  difficult,  it  is  in  this  region  that  the  answers  to  many  of 
the  outstanding  questions  lie. 

Information  is  required  on  the  dynamic  bt  '>ovlor  of  earth  materials  (rocks  and  soils)  over 
all  ranges  of  stress  from  the  vaporization  pressure  to  the  region  of  infinitesimal  strains. 

These  data  must  be  collected  for  a  wide  variety  of  materials.  The  search  for  correlations 
between  these  dynamic  properties  and  convenient  laboratory  tests  should  be  carried  on  vigorously. 

In  the  analysis  of  data  relating  amplitude  to  yield,  the  effect  of  the  properties  of  the  medium 
and  the  explosive  material  cannot  be  Ignored.  This  is  particularly  true  lor  chemical  explosives. 

In  addition.  It  is  only  recently  that  the  Influence  of  ambient  stress  level  or,  equivalently,  the 
depth  of  burial  In  a  particular  test  site,  has  been  included  in  a  quantitative  manner. 

Because  the  relation  between  amplitude  and  yield  is  diagnostic  of  processes  at  the  source, 
detailed  measurements  should  be  made  at  every  opportunity.  The  uncontrollable  variability  of 
parameters  affecting  amplitude  will  make  the  result  of  any  single  experiment  of  doubtful  value 
for  wide  application.  The  results  of  many  experiments  may  reveal  the  correct  formula.  It 
seems  that  problems  of  Instrument  placement  and  departures  from  homogeneity  of  the  earth 
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between  the  source  and  the  Instrument  may  be  as  much  a  source  of  variability  as  the  conaltions 
at  the  source. 

The  dependence  of  the  amplitude-yield  relation  on  the  part  of  the  specirum  *hnt  na  '  been 
observed  on  a  particular  selsmojr-.am  must  not  be  overlooked.  Apparent  scale  effects  In  sroing 
from  small  to  large  yields  may  be  explained  In  terms  of  the  fi  r'^.ueucy  at  which  the  spectrum 
reaches  Its  peak  value  and  the  passband  of  the  surrounding  earth. 
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Appendix 

PRINCIPLES  OF  SCALING  AS  APPLIED  TO  UNDERGROUND  EXPLOSIONS 

1/3 

Throughout  this  report,  the  us«  of  Ihc  cube  root  of  tho  yield  W  ■  lao  *"000  to  fix  ihe 
scale  of  the  explosion.  A  discussion  of  the  principles  underlying  scaling  procedures  will 
be  presented.  Several  authors  have  treated  modeling  and  t'»e  theory  of  similitude  thoroughly, 
but  two  In  particular  have  devoted  themselves  specifically  to  Iht  ground- shock  problem:  Lamp- 
son  [30]  and  Parkin  [81].  Parkin's  paper  contains  an  extensive  bibliography. 

Parkin’s  discussion  is  more  comprehensive  because  he  considers  scaling  from  one  medium 
to  another  as  well  as  scaling  from  one  set  ol  dimensions  to  another  in  the  same  medium.  He 
also  Includes  scaling  for  visco-elastic  and  elastic-plastic  substances,  in  addition  to  elastic 
media.  Only  the  case  in  which  a  single  medium  is  Involved  will  t'e  Included  here.  From  another 
viewpoint,  it  is  assumed  that  in  going  from  one  explosion  (treated  as  a  model  of  all  explosions 
in  a  given  medium)  to  another,  the  velocity  of  elastic  wave  propagation  does  not  depend  on  the 
scale  of  the  experiment. 

The  scale  factor  is  designated  by  S,  and  relates  lengths  in  two  experiments  by  >  SLj. 

On  the  assumption  that  the  velocity  of  propagation  is  Independent  of  the  scale,  the  time  will 

scale  in  the  same  manner;  l.e.,  Tj  =  STj.  The  third  dimension  needed  to  describe  the  quantities 

involved  in  the  analysts,  mass,  is  treated  by  using  density  as  a  parameter  characterizing  the 

3  3  3 

medium,  and  then  M2  *  pS  L  ,  or  M2  =  S  M^^. 

These  fundamental  dimensions  can  be  Introduced  into  a  dimensional  analysis  of  any  of  the 
quantities  of  interest  In  order  to  determine  the  ccaling  relations.  For  example,  stress,  with 
dimensionB  ML'^T'^,  scales  to  S^M  s'*  L"^  c  ^  L’^  *  ML’^T*”^,  so  that  at  the  scaled  distance 
SL,  at  time  ST,  the  same  stress  is  predicted  as  was  observed  at  L  and  T  in  the  original  system. 
A  similar  result  follows  for  particle  velocity.  On  the  other  hand,  particle  displacement  u  at  a 
distance  r  would  appear  as  a  displacement  Su  at  a  distance  Sr  in  the  new  system.  It  must  be 
understood  that  this  analysis  tells  nothing  about  how  the  quantities  vary  with  distance;  one  does 
not  find  the  value  at  a  distance  other  than  the  scaled  distance  from  such  considerations. 

The  choice  of  the  scale  factor  is  arbitrary.  In  the  case  of  an  underground  explosion,  the 
scale  of  the  experiment  is  certainty  set  by  the  yield.  For  a  spherical  charge  of  weight  W,  of 
constant  density,  the  radius  of  the  explosive  charge  becomes  a  length  representative  of  the 
scale.  For  this  raason,  a  length  numerically  equal  to  the  cube  root  of  yield  has  been  selected 
as  the  scale  length.  If  charge  weight  is  thought  of  as  equivalent  to  energy,  tho  cube  root  of 
the  weight  obviously  does  not  have  dimensions  of  length.  Actually  this  is  no  problem  as  long 
as  the  origin  and  function  of  this  numl>er  are  understood.  However,  in  order  to  remove  this 
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apparent  Inconsistency,  Parkin  [27,  81]  has  sujtpested  dividing  the  explosive  energy  W  by  a 

3 

quantity  having  dimensions  of  stress.  A  quantity  with  dimensions  of  L  is  obtained.  He  uses 
the  value  of  Young's  modulus  as  a  characteristic  quahtity  for  an  elastic  me  Jium  with  the  dimen¬ 
sions  of  stress. 

Table  11  Is  taken  from  Lampson  [30]  and  summarizes  the  manner  In  which  frequently  oc¬ 
curring  quantities  scale.  Parkin  gives  a  more  general  table,  valid  for  going  from  one  medium 
to  another  as  long  as  Poisson's  ratio  Is  the  same  In  both. 

TABLE  11.  MODEL-LAW  RELATIONS 

The  ratio  of  scales  Is  S  which  represents  a  multiplication  of  all  linear  dimensions  of  the 
experiment  by  this  factor,  where  S  - 

Quantities 
Quantity  Comparable 
Dimensions  In  at  a 

In  Scale  New  Constant 

Quantity  Symbol  Dimensions  New  System  Factor  System  Value  of  r/S 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Length 

L 

L 

SL 

S 

SL 

L/S 

Mass 

M 

pL^ 

pS^L® 

s' 

s'm 

M/s' 

Time 

T 

T 

ST 

s 

ST 

T/S 

Force 

F 

MLT"^ 

S^M-^ 

s' 

s'f 

F/s' 

Energy 

E 

ml^t'^ 

S^ML^T'^ 

s' 

s'e 

E/s' 

Pressure 

P 

ML'^T"^ 

1 

.« 

p 

P 

Velocity 

V 

LT"* 

•  LT-l 

h '  y 

V 

V 

Total 

Impulse 

I' 

MLT"^ 

s^mlt"^ 

s' 

s'r 

I/s' 

Impulse  per 

Unit  Area 

1 

ML'^T** 

sml'^t’^ 

s 

SI 

I/S 

Displacement 

L 

SL 

s 

SD 

D/s 

Acceleration 

a 

lt'2 

S^^LT*^ 

s*' 

S’^a 

Sii 
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California  Matlta'a  of  TaaRoalogy.  •aUmaloglcai  uhoraiory.  ISO 
North  Ran  Rahal  Avanua,  Faaadana,  Caliromla 

ATTNt  Dr.  NafoBanioff 

California  Snotlmia  of  Taahnalagy.  Balamaloglaal  Labaraiory.  ISO 
North  San  Raraal  Aranva,  Faaadana.  California 
ATTNi  Dr.  1.  Knepafr 

CallfamlB  Inatlnita  of  Taahnolacp.  lalainalagtaal  Labaraiory.  tso 
.North  San  Rataal  Pv*"**#,  Fieedana,  C*!!f‘?r>t« 

ATTNi  Dr,  Frank  Fraaa 

Vhlrorally  of  Cailfarnla,  looranaa  Radiation  Labaralary.  tJvarmara. 
Califamla 

ATTNi  Dr,  Roland  F.  Narbal,  Dr.  Charlaa  B.  Vlalat  (tt 

Vnivaraity  af  Califamla,  Sarippa  Mtlltata  af  Oratnagrlphyi  to  Jolla. 
California 

ATTNt  Dr.  Wallar  Manh 

Vnhitrally  af  California,  Salatnafraphla  tNallan,  Barttaloy  4,  Calthir* 


ATTNi  Frafraaar  F.  Byarly 


Camofto  fottflvia  af  WaahJnften,  DaparWK  •*ioatsidl  Magno- 
tlom,  SMI  rf^g  Braneh  Road,  N.  W.,  Waahlngton,  U.  C. 

ATT'.  J.  B,  .Siabihart 

Caotary  Caaphyaleal  CorparallM.  F.  O.  Bon  C.  Admiral  Stallen. 
Tolaa  IS.  Oklahoma 
ATTN:  R.  A.  Brodp^ 

Caatary  Gaapbyolaal  ^’‘oaaratlan,  F.  O.  Boa  C.  Admiral  Station, 
Tataa  IS,  Alahama 

ATTN:  J.  e.  A.  T.  Vaadbom 

V.  t.  Canal  and  Onadaile  Sarray.  Camp  Mareary,  Navada 
ATTN;  T.  &  Floraa 

Daparmeaat  af  Cammaraa,  Caaai  aad  Oaadatle  Santay,  Oaephyataa 
Dhrlalea,  VaaRb^ian  tt,  O.  C. 

Al*^*'-  Or.  Daaa  t.  Cardar 

Dapai  mill  ad  Cammarn.  0.  i.  Canal  and  Oaadatle  flantay,  Oaophy 
oica  Olalalan,  WaMilafian  SI,  D.  C. 

ATTN:  Mr.  J.  iacAM 

Dapartinmi  ad  Cammaraa,  Cm*l  and  OaodaHo  Sarray.  Oaapfiytlca 
DHriaioa,  tPaafclng»n»*  Ik,  p,  r. 

ATTN:  Mr.  Loanard  MarpAy 

CaMlnanml  DU  Campnny,  Fl  'd.  Oraorr  StS,  Ponca  City,  AMhema 
ATTN:  Dr.  Jidm  M.  Crawlfrrd 

Dafmaa  Aiamle  tiBPon  Aganoy,  Flald  Cammano  Comma.idaf.  %ndla 
Baaov  ABapaatNoa.  No*  Maateo 

DafOnaa  Alomie  Sapport  Aganay.  Baadgaariara,  VaaMagian  tf,  D.  C. 

Arm  M»o.  uvip 

Draaaar  DoetronJea,  MB  Dhrlolaa,  F.  O.  Boa  ISIIT,  lloaoten  ST. 

Taaaa 

ATTN:  Mra.  MUdrad  Baami'a."  UbrarUn  (t| 


Bdgarion,  Oarmaahaaaan  and  Orlar,  Dlraeiar  af  Raamrah,  lid  Broeh* 
lino  Araviaa,  Boa  lew,  MaameWaatta 

ATTN:  Clyda  DabbI# 

Bdtanea,  Oarwiibaaam  and  oriar.  fceorporaiad,  ISO  Broehllna  Atra*' 
noo.  Bearnn,  Maaoachoatita 
ATm  F,  T.  Birabela  *** 

Blaclr*  v.*ghanlea  Compoay,’  P,  O.  Bon  srs.  Aaatin  M,  Taaot 
ATTNi  Df.  Frad  N.  Merrla^V^ 

jiT*.*' 

BnglMaar  Baaaareb  aad  Oioi wt  tobaratary,  Minor  Ataailan 
Branch,  Far!  Balvoir,  Vlrglnm  - 
ATTNi  Mr.  B.  B.  Dwamli 

P. 

Bnglnaarlng'Fhyalea  ComaaiBN'^^^'^  Randolph  Mraal,  RaekvlUo, 
Maryland 

ATTN:  Mr.  VMeant  CaabJng^)' 

Oanarat  Atranlaa  Carimrailai^|)na  Bala  Aranua.  Bala^Cynvyd, 
Fonnaylmnia 

ATTN:  Mr.  Joaapfi  T.  Vnda^^,  tn 

Qanamidaaphyalaal  Campan^TSO  Nouaian  Club  Building,  Houetan, 
Te*aa  , 

ATTNi  Dr.  I.awla  M.  Mall-^>t*' 

Qaaphyalaal  and  Polar  Raaaarc^-^**'**^  Highland  Road, 

Madiaan  S,  WlaennaM 

ATTNi  Frafaaaar  Oaarga  Baallard,  Dlraeiar 

Oeaphyaleal  Saralea,  Inaarperalad,  F.  O.  Baa  SI0S4,  Allaa  SI, 

y 

*TTMr  Dr,  0«  MiMar, 
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The  University  of  Michigan 


OiOtelwtegl  Cir|Wftlan.  WMSO,  |^«  O.  Bm  Mtt,  Cik*lM«a 

«» 

e^rf»ftim,  f.  O.  telB*  H.  T«iM 

Am<i  Or.  W.  U.  (M» 

OiWmi*  RMMrrh  Cmit,  P.  O.  9m  MVt«  OiliM  t.  Tmu 
ATTTti  Dr.  iMfi  Oiftour  (t) 


HlMlMlppI  iiiMtrtel  Md  TrrftJitHocIWf  RtMtrrN  TWnml— o. 
Dm  ItfT.  Jm*  M  t, 

Arm  Mr.  e.  / .  mum 

mtmrnl  CwHpMiy,  mi  9mM  ^alr  p»m» 


0«tf  P»marm  mi  Drm 


H  OfrlalM.  OMyhy«lAl  Rforeli  Dl> 
trialan.  P.  O.  Dr:«v«r  MSt.  PinaMirgli  M.  FannrytvMlft 
Arm  Or.  T.  J.  0« Quill,  Olr««««r 

tMitrcralty  af  Hawaii  haiHata  af  Oaopl^aiea.  Hwaatala  M. 

Hawaii 

Ami;  Or.  Daa>Caa.  Oipaif  aal  fcrti  8clwwai 
tMraraHy  af  Oltoala,  Oaparfi  mt  al  Mtali^  UtAm.  tIlfciaCa 
ATTH;  Or.  A.  I 


Arm  i.  r.  i 


mmrni  mtmm  rmiil  .Ida,  Mil  Caaa*i*ftilea  AvMaa,  H.  «NaR. 

,  0.  C. 


OapartMaai  ai  kfariaa.  r.riftvA  a/  Craaial  ihMtaa.  O.  A  QaaiaalaaJ 
•array.  TItO  Waat  l#ia  Mraai.  lAhtarauti  11.  rr<^r«iaa 
Arm:  Or.  OaarfaKallar  ^ 

Dapanaiaai  af  Maarlar,  Hraaali  af  Crwaial  •»w^.aa.  0.  L  Qaalagiaal 
•array,  ¥•••  Haat  lata  iKrtr,  i^aavcwd  li,  CalaraAt* 

Arnii  UalaC.  Oaklaar 

0apar«MKt  af  fertarlar,  0.  A.  Oaataflaal  larw^f.  Daarai  FaAaral 
Caaiar.  Daarar  tf,  Calaraia 

ATTH:  Or.  H.  RaMli 


Oak  UkaraMrIaa.  LaaMfiaa  t1.  Maaaaakaaaiia 
Arm  Mr.  Oak  art  intaakif 

^araay  Frakaatiaa  Raaaarek  Oaaaaaaiy,  IttS  Hank  Lawfa  Art 
Talaa  M.  OklaNoM 

ATT!ti  Mr.  P.  h 


JalBi  earraO  OMraraHy,  Oaraiwii  II,  QMa 
Arm  Oav.  a  f.  Olitaakaaar,  A.A. 

Jaka  Carratl  PklaaraOy,  Claralaal  II,  Okla 
Amt:  Or.  MmH  Waltar 

Laaiant  Oaalaflaal  Okaarwaiary,  CakiiaklB  Oktraralty,  Oallaalaa, 
Haw  Yarfc 

AtTKi  or.  Maarlea  Oafcif 

Uatani  Gaateftaat  Okaarvatary.  CalaMkla  OaMaraliy,  Qallaalaa, 

Haw  Yarh 
ATTN:  Or.  H.  I. 


■i — -n*  Oaalaflaal  Okaarvaaary.  Calawikto  MilrwraOy,  VatlasAaa,  Ha* 
Tark 

ATTH;  Or.  tmk  B.  Oltrar 

Maaaa*liaaatta  MaiHaia  af  Taakwalafy,  Daparimawl  af  tartk  IH'tMsaa, 
Camkrll^  11,  Maaaa4kaaatta 
ATTM  Orafaaaar  A  M.  Haafaaa 

Malfar,  Maarparaial,  Afollal  •elaaaa  Otvlataa.  11  Oalaa  Kraal, 

Watariawa  rt.  Maaaaakaaatta  .  i. 

ATTM  Mra.  LarraKa  Matara.  Ukrarkwi 

Maliar,  MaaroaraMI.  TaaM  Hal  laMraiatika  Caatar,  looe  ArliMftan  < 
Haalavarl.  raUa  Ckarak,  virflNla 


ATTHt  r.  D.  VK 


.h’ 


OaHaraHy  af  MloaiMi,  Matllata  af  TaakMlaei.  iakaal  at  MKaa  awi 

Maialhirfyi  MKaaaialta  M,  MkwaaaM 

ATTNi  Oralaiaaf  tOraM  M.  Maaaay,  Oaafkyataa 


>,  Offkia  af  Haral  Oaaaarik,  Waakkifiaa  11,  D.  C. 


Hkay  BaetraaMa  lakaa alary.  U.  A.  OIractar.  9m  Dlaga  U,  Califkr- 

ala 

Aim  Mr.  $.  MaMfitea 

Mary  OaMaMflaat  Palaaa«  takaratary,  Dwratar,  laa  rrr«elaa*  M. 
CalHantK 

Arm  Hvi  t^awa  Caafar 

itoay  AadMHvHa*  DafM  Ukaratary.  Hmiara  Oakit,  OalMkif  Ilf, 
laa  riaaalaaa,  Cailtarala 
Arm  Or.  1.  Oaaankaiaa 

Mary  OaMalaflaal  Dafkaaa  la  karatary,  04f  aatar,  Im  fraaaUai  M, 

CallMfwla 

Arm  Mr.  HkiMmHO 

Cr^  Mata  Callafa.  Oafariaiaat  af  OiaaaapiiMi,  Oarvallla.  Oragaa 
Arm  Dr.  9m  Oarf 

Orafaa  Mata  CaUaga,  DafaTtaaaal  af  Oaaaaafrayky.  CarratUa,  Orafaw 
Arm  Dr.  P.  DakUatar 

Tka  Bklrarally  af  OkOkaaia,  Oaaaarak  baMHata,  Haraiaa,  CkOkf  a 
Arm  Dr.  Harawa  OWiar 

PaaaayOaalB  Kata  Dkh^alty.  Payartaiaat  af  Oaaykyalaa  aal  Oaa* 
akaailalry,  ttO  Mlaaral  lelMiera  OalUkif.  Ikitewrally  Park,  Paaaayl* 
awOa 

Arm  Dr.  D..  r.  trnmU,  Jr. 

RaMaCarparatlaaaf  Awarlea,  Caatlaa  I,  Havfaraay 

Arm  Mr.  O.  w.  K.  Kkif 

HaMa  Carparaftaa  af  Aa^laa,  DarM  SaraaO  Oaaaarak  Caatar, 
Prlaaataa,  HarJaraay"* 

Arm  Dr,  D..  A  MaCay 

Oaal  CaryaratHa.  tka,  ITOO  Mkki  Kraal.  Mala  Maalaa.  CaMfkrala 

Atm  Dr.  RIaharl  Uttar 

Oaaaaalaar  Palyiaatatta  Katitvta,  Tray,  NavYatk 

Aim  Dr.  OMaMlKaiav 

Oalraralty  af  OUla  totoM.  Kkiyaiaa.  Ohala  lalaal 
Amt:  Prafaaaar  Kaaa;.<^'' 

.  0  Ha  Oalraralty,  Haaataa  1,  Taaaa 

*Arm  Prafaaaar  I,  Ct  DaOramaaekar,  OayarKta*.*  af  OaalaB 

Oaiaa  Air  Daratoytnaai  Caatar.  MaatMaarlara.  Orlfttaa  APt,  Ma« 

Tatit 

Aim  lit  u.  9mm  N.  DPitamlafar,  tiMP  ID 

Mafia  Caryaritlaa.  OWItlaa  Tono.  MafU  Oaaa,  Alkuauar^aa,  Na* 

Maalaa 

ATm  Mr.  0.  A,  Pawlar 
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Institiit*  of  Science  and  Technology 


Iftndto  C»ri>or«tlm,  Otvlstan  SI1S<  •ftndl*  ■•••. 

MtKlee 

ATTNt  Mr.  M.  I.  M*rrm 

iMidia  CerporBllM,  DIvlalM  TSll.  SaAc'U  AlMumi^w .  Nm 
M«xleo 

•mMM  Cai>«ratlM,  iMdta  MiM,  hnmumrtm.  Him  MmIm 
ATTM;  WUltaiil  R.  P«rr«t 

tMMM  C9rpofi\m,  OHrlifoa  fttl.  teiMM  Ram,  Aftm«*rRM.  Mw 
MmIm 


•mU  OmlopmAiM  rumfamf,  R.  a  Bm  Ml,  Hmmmi  1.  T«n» 

Amt:  Or.  K««rfiuMi 

MrtHMlM  lAilv^rilty,  OtnilBi— i  M  OMtoflp.  teOM.  TAM* 
ATT?f;  Dr.  OxfiM  T.  Ilirria 

•pMA-OAMral  CarRarattM.  tTf  f  laa^r  llrati.  !,  A»t**aniia 

Arms  Mr.  A.  P.  AlWaaM 

%>ai  OawararOarparatiaa,  tfV  RMaar  Rtrax,  0l«»ar'4»  t,  Ca*<tarala 
Arms  Or.  OlMV  U  »*«B 

Ipaaa  qimiral  CartaratiM,  tff  Riawar  tiraat,  Olaaiala  1,  CalOarala 
Amts  Dr.  R.  R.  Otiaaa 

ItaafarA  RaaaaraA  Madtuta*  Maala  Rarti,  CatUamla 
Arms  C.  M.  AMa* 

twafirR  RaiMrah  feMMltaM.  Maala  RaiR,  Oatttdrala 
Arm.  Or.  DarU  DaniataM 
SiaMara  Raaaarah  Mttliata.  Maala  Rait^  Calliarata 
Arms  Or.  Oaarfa  t.  OaraU  ^ 

•MtaM  R.Mrtk  kiniMi.,  ewdur  UkwkiM'IM,  MMI»  M«t;  OlSl* 
larala  C'r  ^ 

Arms  Or.  lata  O.  tHiaiM 

Riaaf«r4  RaaaaraA  MMiaata.  Maala  CAlMirala 

ATmi  Dr.  Rlatar4M.  Raaaa 

•tatfara  RaaaaraA  Mataaia.  MaaM  Rart.  CiUlarala 
Atms  f.  R.  OrMa 

■iMrfara  Raaaarali  laatHiAa.  Raia  ANa.  Caltfarala 
ATTRl  nr.  AUaaRafaraaa 

■laafara  RaaaaraA  MMUata.  RaaRar  Utaralartaa.  Mania  tarts.  Call- 
(arfia 

ATTRs  Dr.  Ttaaiaa  C.  Raatar.  DIraaiar 
Maalara  RaaaaraA  MaWirta,  Maala  Rarti.  Caltfarala 

Arms  U  M.  taio 

MaafarA  RaaaaraA  MMItaa.  DaUiMf  IM,  Maala  tarti.  Caltfomia 

ATTfti  Dr,  Rakari  D.  Valla,  Ir. .  Dlraatar.  RRjralaa  fMvtaton 

Maafard  Unkaraliy,  NliR  taarQr  tayalaa  tataraiary.  fltanfara.  Call- 
famU 


91.  Laala  tialraral'  .  Tlia  laMitui*  af  *l8cla».i:s;  Tt*>  nnbvi '•.f?!, 
•t  taala  1.  Mlasaari 

ATTNs  Or.  R.  R.  Ilaifirleh 

M.  Laala  OahraraHv.  nr>  -tei«e.;*>f  Oaophyalea.  Stl)  Oliva  Utraat. ' 
M.  Laaia  •«  Miaaar-. 

ATTiai  Dr.  Vu.lli*  atkaAi^.  S.J, 


•*.  I*.  O.  ^  Daisaa  9S,  Taxaa 


ATTK:  R.  A.  Araatt  (f) 


ATTRi  Dr.  WatffiaflM.  N.  taaaNky 


Rl.  laala  iMIvaraHy,  Tta  MallHiia  al  Taataalagif.  Itit  OHva  Rtraai. 
as.  lauia  I.  Mlaaaarl 


ATmi  Of.  CarlRIaalMfar 


Oairavalty  af  Tanaa.  Aaatia  If.  True 

ATT1R  RiMaaaar  W.  T.  MaalUbr^f*^*  Dapartmant  at  Orelofjr 

Drrtaa  tar«l  telaaeaa  Dlvlataa.  R.  O.  Baa  314.  Alanndrto.  Vlrfbito 
Arm  Mr.  P.  B.  Cakar  n» 

DaNaA  BaetraDyaaailea,  taaarpo»>ta4.  100  AllanAala  hwaJ,  taariOana 
CalMarair 

Amts  Mra.  Daratky  A.  Allan  U1 

MoHr#  Sia«*raDrtiaMka.  laeMfMV  tOo  Atiaatato  RaaO,  taaadana. 
Calllatala 

Amii  C.  B.  Rartaa 

OaNaO  BaetraOjnMlaa.  Maaraaratatf.  too  AUaniala  Raad.  taaaAaaa. 
CalDamla 

ATT3R  Mr.  Allaa  M.  Ra«.  Ir. 

D.  i.  OaaMflaal  tarvay.  Ml  MIMlafsaM  Raad.  Mlnia  tarts.  Callftor* 

ala 

ATT1R  Utaarlaa 

Oalraral^  at  Dials.  Oapartaiaat  at  Oaaphyatea,  Mlaaa  BatMtaf;  MH 

Uisa  City,  ttah 

Arms  Rrataaaar  Kaa  Carts 

■  ■  .V*-' 

RfooMi  Otalrsrtlari,  Waliaa  tl.  MalaMlkaOatta 
At^  italraaifiaalal  Unataa.  tii. 

Wioaaaala.  OaapRyalal  al^  iNlar  Raaaarek  Csmtar. 
tMA4  tMtrtraliy  Avansia.  MaAtaan  I.  sliiaaaaifs 

lUrlta  DmataAliy.  CMclaaatl  T.  OktO 
ROtariaA  B.  A.  Drtalay.  8.1. 

a^atale  W  ClaaAa  Bataairt.  taria  Rranaa 
A\*rNi  Rrataiaar  T.  RoaarO 

Aratle  laailtata  at  itarDs  Amarlea,  3411  Radpails  ftraat,  Maniraal  tl. 
R.D**  CiaaOa 

Arm  Mr.  lata  C.  RaaA.  taaavttva  Dlraatar 

laaa  Latal.  Alpaaa.  10  Raa  ClaaO#  BamarO.  tarla  V.  Rraaaa 

tatrtaalty  at  taaiiaiataaan.  ■aaliaieen,  R.aliai«ha<«an,  Canaia 

Arm  Dr.  lamaa  Mataalay 

fal^MMlatlaal  tataratary,  Oppaalki  taaOan 

Arm  Dr.  Martsaa  Bitls 

VisIvaralD  af  Taraaia.  Daya^taant  at  RHyalea.  Taranta  I,  CaiiaAa 
ATT^fi  Rrataaaarl,  T.  Wllaen  111 

WakmaM  Mailtala  at  Atiaaaa,  R.  o.  Baa  ll,  Ralsavasis,  taraal 
ATTNi  Rralaaaar  C.  L  tatsarla 

Ualvarllty  at  Wltatsarartsrt.  fesailtuia  at  Rarfsar4  Rrkt.  Qaafhyatsval 
Raaaarali.  lotanistakuri.  Aassih  Africa 
ATTHi  Dr.  A.  U  Milaa 
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Th#  University  of  Michigan 


Advanced  Research  Projects  Agency 
The  Pentagon,  Washington  3S,  D.  C.  (-S^ 

Advanced  Research  Projects  Agencf 
The  Pentagon,  Washington  29,  D.  C. 

ATTN:  Mr.  RuaseU  Batrd 

Advanced  Research  Projects  Agency 
The  Pentagon,  Washing^  29,  O.  C. 

ATTN:  Mr.  R.  Black 

Advanced  Research  Projects  Agency 
The  Pentagon,  Wash'nfton  29.  f).  C. 

ATTN:  Mr.  Donald  CUiiienU 

Advtneed  Research  Projects  Agency 
The  Pentagon,  wascir^on  ye  ;i.  C. 

ATTN:  Mr.  MaJ.  R.  Rarrls 

Advanced  Research  Projects  Agency 
Tlw  Pentagon,  WashlngtM  29.  D  r. 

ATTN:  Dr.  Jack  Rulna 
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